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PREFACE 


The investigation documented in this report constitutes part of the lunar 
roving vehicle (I.RV) research conducted by the Advanced Lunar Studies Team 
at the Jet Propulsion Laboratory. This study was performed to develop and 
provide a better understanding of mobility concepts on soft sloping terrains as 
applied to LRVs, 
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ABSTRACT 


A general solution is given to the mobility performance problem of a 
power- driven rigid cylindrical roller climbing a semi- infinite soft soil slope 
with uniform velocity. The roller axle is subjected lo vertical and pull force 
components. A gravitating, cohesive-frictional soil is considered. Its appli- 
cation to lunar and planetary locomotion is emphasized. The mechanic s of soil- 
roller interaction is described and solved, considering stresses and velocities, 
as a mixed boundary value problem. Kotter's quasi- static equilibrium equa- 
tions arc connected to a plastic stress configuration satisfying Shield's velocity 
conditions along the characteristic lines. Solutions of the equilibrium equa- 
tions yield the driving torque, slip, sinkage, and soil- roller interface stresses. 
Driving power requirements and thrust efficiency are determined. 

A general concept of safety factor against immobilization is introduced. 

A computer program for the soil-wheel interaction performance |SWIP) was 
developed and limited applications of this theory to rigid wheel tests on hori- 
zontal terrains indicate very reasonable agreement. The method is also 
applied to the Apollo and Lunokhod- 1 lunar roving vehicle wheels. Pari I, as 
published, presents the basic and necessary conditions satisfying the limiting 
equilibrium and velocity equations. Part II, to be published separately, will 
provide the concepts of sufficiency asserting the completeness of a given solu- 
tion and the computer program. 
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T. SUMMARY 


The use of wheeled roving vehicles in future 
lunar and planetary explorations will require con- 
sideration of two technical problems: (1) the ex- 

pense and difficulty of performing lunar soil sur- 
face tests to derive, experimental coefficients 
which define soil-wheel performance, and (3) the 
lack of the basic fundamentals and experimental 
background to predict either on a theoretical or 
empirical basis the mechanics of lightly loaded 
rolling devices on soft horizontal or sloping 
terrains. 

This study presents a method of solution to 
the two - dimensional problem of a power- driven 
long and rigid cylindrical roller moving up a gen- 
erally sloping soil surface with uniform velocity. 
The roller axle is subjected to the combined action 
of a driving torque, a vertical load, and a pull 
force parallel to the terrain slope. All practical 
ranges of soil mechanical properties are con- 
sidered for either lunar or earth-based soils, in- 
cluding soil friction, cohesion, and gravitational 
effects . 

In principle, the soil-roller interaction study, 
as would one for a wheel, points to two basic as- 
pects. One corresponds to the operational condi- 
tions of a roller whose mobility is always guaran- 
teed, particularly for low contact pressures and 
reduced sinkage. In this ease, the design objective 
for mobility is mainly concerned with an efficiency 
optimization problem. The other aspect relates 
to the particular limiting state in which immobili- 
zation occurs when the soil thrust capacity has 
been reached due to the combined action of soil 
weight and applied roller loads, in this case, a 
margin of safety against immobilization, rather 
than efficiency, constitutes the constraint design 
factor. 

This study indicates that the soil- roller solu- 
tion represents, to a reasonable degree, a 


satisfactory basis for estimating the performance 
of rigid wheels under .similar driving conditions. 
This is particularly true when the soil-wheel con- 
tact pressure level is low and the soil failure pat- 
tern develops mainly in the fore- aft direction 
rather than in a lateral mode. This characteristic 
is expected to occur predominantly for vehicle 
locomotion on the lunar surface, as evidenced 
from a photograph of wheel tracks left by the 
Soviet Lunokhod- l (Fig. 3 of Ref. 1). Obviously, 
the soil-roller solution does not entail the solu- 
tion of the wheel problem 'per se'" since the wheel 
is a finite representation of a roller. However, 
the roller approach is a necessary step which will 
permit the elucidation and disclosure of some im- 
portant aspects applicable to the soil-wheel inter- 
action problem. In particular, these aspects 
refer to : 

(!) The mechanics of slope climbing under 

various types of soil materials, and load- 
ing conditions of self- p rtipuls ion or pull 
forces. 

The role of slip as a kinematic factor 
affecting the wheel mechanical perfor- 
mance and locomotion efficiency. 

(3) The basic principles concerning how' soil 
thrust is generated and motion sustained 
as a continuous mechanical process. 

If all these factors can be derived for a roller, 
the conclusions may be directly extended to wheels 
or otherwise approximated by appropriate model- 
ing techniques. In this context, the roller analogy 
is more akin to the soil-wheel interaction problem 
than to using flat plate tests since in the former 
the mechanics of rolling is reflected at all slip 
levels as it. actually occurs for a wheel. 
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The solution approach consists of selecting a 
soil-roller failure configuration in accordance 
with experimental evidence on horizontal terrains. 
This failure pattern is generalized to sloping sur- 
faces. A compatible velocity field is defined that 
satisfies all velocity requirements. Then the 
governing stress equilibrium equations are solved 
along the velocity characteristics (sliplines) in 
connection with the remnant stress boundary con- 
ditions. The method of solution is implemented 
by a computer program which evaluates the basic 


soil-roller performance parameters as given by 
the torque, slip, sinkage, and interface stress. 
The concept of specific rolling energy dissipation 
is generalized to slopes, and the driving power 
requirement is determined under quite general 
conditions of load combinations and terrain slope. 
A general mobility safety factor definition is in- 
troduced and calculated for a soil-roller system. 
This definition applies also to wheels and covers 
the whole spectrum of rolling with slip up to and 
including immobilization. 
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II. BACKGROUND - LUNAR MOBILITY 


The current status of soil-wheel interaction 
as a design tool for planetary vehicle exploration 
is rather unsuitable due to the lack of pertinent 
quantitative tests and design information that is 
normally available for on- earth vehicles. New 
theoretical concepts are required to predict vehicle 
performance requirements, especially with regard 
to safety and power needs when traversing lunar 
soft terrain and steep slopes. In this context, a 
pertinent survey of the state-of-the-art on off- 
the-road vehicle design was made and has led to 
the following conclusions: 

(1) The present status of the theory of soil- 
wheel interaction applies exclusively to 
horizontal or gently sloping terrains. 

Wheel design parameters and vehicle 
performance are definitely connected 
with the soil mechanical properties and 
the nature of lunar topography. Con- 
sequently, lunar vehicle slope climbing 
and traversing capabilities are bound to 
constitute a controlling design factor in 
terms of mobility safety and locomotion 
energy requirements. 

(2) Current on- earth vehicle mobility experi- 
ence refers mainly to high soil-wheel 
contact pressures and dictates that heavily 
loaded vehicles cannot practically negoti- 
ate soft slopes higher than approximately 
25 deg. Instead, lunar roving vehicles 
will operate under comparatively low 

soil contact pressures (approximately 
6. 9 X 1()3 N/m^, 1 psi) and may require 
climbing slopes steeper than 25 deg. It 
is not known, using present mobility con- 
cepts, if a vehicle can safely operate 
under conditions of low pressure and re- 
duced gravity on lunar crater slopes near 
limiting equilibrium where excessive 
sinkage and loss of soil stability support 
may be imminent. 

(3) Present approaches to derive basic per- 
formance parameters of soil-wheel inter- 
action resort to separate sinkage and 
horizontal shear deformation plate tests 


performed on potential terrains of vehicle 
operation. On this basis, a comprehen- 
sive semi- empirical theory for off-the- 
road mobility was developed by Bekker 
(Ref. 2). The inherent uncertainty of 
this approach relates to the fact that the 
character and distribution of soil-wheel 
interface stresses derived from flat plate 
tests cannot be reliably translated into a 
roll on a sloping surface. 

(4) There is insufficient evidence to what ex- 
tent the combined or independent action 
of friction, cohesion, and gravitation in- 
fluences soil-rolling thrust performance. 
Experiments and analysis by Costes, et 
al. (Ref. 3), indicate that the lunar soil 
surface mechanical properties relate to 
a frictional- cohesive behavior, 

The use of wheeled roving vehicles in future 
lunar and planetary explorations will require con- 
sideration of two technical problems: (1) the 
expense and difficulty of performing lunar soil 
surface tests to derive experimental coefficients 
which define soil-wheel performance, and (2) the 
lack of the basic fundamentals and experimental 
background to predict either on a theoretical or 
empirical basis the mechanics of lightly loaded 
rolling devices on soft horizontal or sloping 
terrains . 

The broad objective of this study is to estab- 
lish a basis for a consistent and unified approach 
to the soil-wheel interaction problem to estimate 
the performance of lunar roving vehicles travers- 
ing generally sloping soft surfaces. To this end, 
basic soil mechanical concepts, applicable to the 
lunar environment, are incorporated in analytical 
expressions to derive torque, energy, interface 
stresses, operating slip, and sinkage values as 
may be applicable to lunar mobility requirements. 

These analytical expressions can then be com- 
pared with the results of controlled experiments 
either on earth (Ref. 4) or on the moon (Refs. 1 
and 5), using, for example, the Soviet Lunokhod- 1 
lunar roving vehicle mobility system. 


JPL Technical Memorandum 33-477 


3 



III. SOLUTION APPROACH 


The solution approach to the soil- roller inter- 
action problem is based on the application of the 
theory of plasticity to soil mechanical problems. 
The governing differential equations of equilibrium 
(Ref. 6) and velocity compatibility along the slip- 
lines (Ref. 7) are solved as a mixed boundary 
value problem connecting stresses and velocities. 
This problem is concerned with a quasi- static, 
steady- state, two-dimensional plastic flow with 
soil assumed to behave as a rigid perfectly plastic 
material. The soil is granular in character with 
strength properties defined hy the Mohr- Coulomb 
theory of failure, which depends on the soil angle 
of internal friction cj> and cohesion c. The soil 
strength properties are isotropic and homogeneous 
throughout the plastic domain up to and including 
the soil- roller interface. ^ Soil self-weight is 
specifically considered. 

The problem of rolling contact between a rigid 
towed cylinder on a plastic deforming horizontal 
half space was investigated theoretically by 
Marshall (Ref. 8) for a Tresca material, applying 
a perturbation method of solution. Dagan and 
Tulin (Ref. 9) studied the steady flow of a rigid 
plastic clay beneath a driven cylindrical roller on 
a horizontal half space, using the method of slip- 
lines. Experimentally, Boucherie (Ref. 10) ob- 
tained photographic flow patterns produced by 
towed and driven rigid wheels on a supporting bed 
of packed metal rolls (Fig. 1). Wong and Reece 
(Ref. II) performed a series of tests on level 
sand surfaces, acted upon by rigid rollers and 
wheels for various loads, slip, and skid combina- 
tions. Using photographic techniques, they pre- 
sented and commented in detail on the character 
and patterns associated with soil failure due to 
rolling loads. 

To date, all known complete solutions of un- 
contained rigid perfectly plastic flow problems 
consist of selecting an appropriate slipline stress 
field and then verifying whether or not the bound- 
ary and the field velocities are also satisfied 
throughout. If not, a new stress slipline pattern 
has to be tried (Ref. 12). Obviously, this pro- 
cedure is rather limited regarding the boundary 
velocity functions which can be considered other 
than simple uniform ones. 

Here, instead, an inverse procedure has been 
adopted. First, a compatible soil-roller velocity 
field of characteristic lines is generated which 
satisfies the roller boundary velocity conditions, 
then to this field is adjoined a stress domain that 
satisfies all the remnant boundary conditions. 

That this is possible is based on the fact that 
velocities and stress characteristic lines are 


coincident (Ref. 7) if the former derive from the 
theory of plastic potential (Ref. 13). To the 
author's knowledge, this is the first time a prob- 
lem of mixed boundary values in the theory of 
plasticity is solved for both stresses and veloci- 
ties by satisfying first the velocities rather than 
the stress conditions. 

This method considerably facilitates the 
search for a complete solution. Also, this ap- 
proach may be found useful in the study of metal 
forming and tooling operations and other soil me- 
chanical problems which are accompanied by vari- 
able interface friction and complex moving plastic 
boundaries. There exist infinite stress field pat- 
terns which can satisfy equilibrium. In addition, 
completeness of a solution is defined only when 
(1) the kinematic compatibility conditions are also 
satisfied, (2) the dilation rate is positive through- 
out the plastic domain, and (3) at no point of the 
rigid domain do the stresses exceed yield. It is 
reasonable to expect that if the problem is well 
posed and a solution is found, the results obtained 
will coincide with a possible matching experiment. 

The sequence of the study reported here is as 
follows: In Section IV, a velocity flow field pattern 
is selected and validated in accordance with experi- 
mental evidence of rolling tests on horizontal soil 
surfaces. These results are generalized to slopes. 
The kinematics of slip and velocity boundary con- 
ditions are analyzed as they relate to the soil- 
roller interface. These velocities are then prop- 
agated throughout the plastic domain. In Section 
V, a summary of the plane strain theory for rigid 
plastic solids is given. The slipline fields de- 
fined in Section. IV are used to determine soil 
limiting stresses. Soil reactions and moments 
are calculated by considering soil weight. The 
quasi- static equilibrium equations are defined and 
the soil-roller interface stresses are evaluated. 

In Section VI, the results of Section V are 
synthesized into a system of equations whose solu- 
tion is the basis for the complete solution of the 
soil-roller problem. The total torque, energy, 
and roller sinkage equations are formulated. A 
general definition of the factor of safety against 
immobilization is introduced. The computer pro- 
gram is applied to soil- rigid wheel test results 
which are numerically compared to illustrate the 
use of the method and to verify its prediction 
capabilities. 

Finally, the theory is applied to the Apollo 
lunar roving vehicle (LRV) and the Soviet 
Lunokhod- 1 wheels operating on the lunar surface 
to estimate their mobility performance. 


*If the soil-roller interface constitutes a material strength discontinuity, particularly if the interface 
friction and/or adhesion properties are different than the corresponding soil-properties, then the local 
orientation of the sliplines is dictated by the relative soil- roller interface mechanical properties. 

There is no theoretical difficulty in accounting for this fact in the analysis in terms of either one or 
both factors, adhesion or friction, as long as the strength discontinuity is represented by a Mohr- 
Couloxnb relationship. 
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XV. SOIL-ROLLER KINEMATICS 


A. Problem Statement 

The roller moves with uniform velocity Vq 
parallel to the original undisturbed soil surface 
(Fig. 2). The surface slope angle ct is measured 
positive in a counterclockwise sense from the posi- 
tive x-axis. The rotational velocity u> (rad/s) is 
positive clockwise. An orthogonal Cartesian coor- 
dinate reference system (x, z) is adopted with ori- 
gin at the roller axle center C. The positive 
z-axis is oriented down, parallel to the local grav- 
ity vector. As a steady- state process, the kine- 
matics of the roller is described by the position of 
its center of instantaneous rotation I(x, z), located 
along a line passing through C and normal to the 
original surface. 

x - - Rs^ sin a (1) 


z = Rs^ cos a (2) 


The translational velocity of the roller center C is 


V c = uRs k (3) 


where s^ ~ slip factor and Rs^ defines the position 
of point I with respect to C. In terms of displace- 
ment, is an index number which, for a full 
wheel turn (2v rad), indicates the axle displace- 
ment as a fraction of the developed wheel perim- 
eter length. For a driven roller, 0 s sh - 1.0. 
Limiting conditions are: 

(1) Pure rotation, = 0 and Vc = 0. 

(2) Pure translation, = 1.0 and Vc = u>R. 

The roller kinematics may also be evaluated, 
based on the knowledge of the rotational velocity 
wand the translation velocity V^, by following the 
definition of slip, 


where V = uR = roller peripheral velocity. It is 
verified from Eq, (4) that, for driven rollers 
s > 0, the roller slips. For towed rollers, 
s < 0 and it is said that the roller skids. Based 
on Eq. (4), s = 100% slip for pure rotation and 
s = 0% slip for pure translation. 

This type of motion description was used by 
Poletayev (Ref. 14) and by Onafeko and Reece 
(Ref. 15) in the study of rigid wheels on level ter- 
rain ( a = 0). Here, this motion concept is gen- 
eralized to slopes where a > 0. The problem is 
now posed as follows: Given a semifinite soil sur- 
face slope a, possessing the unit volume weight Y 
and strength properties defined by the cohesion c 
and the friction angle <j>, loaded by an infinitely long 


rigid roller of radius R, carrying an axial weight 
W and pulling a load P* per unit roller width b 
parallel to the original slope, determine the 

(1) Operational slip factor sj< which defines 
the limiting equilibrium condition that 
permits the roller center C to move with 
uniform velocity Vq (Eq. 3). 

(2) Driving torque M capable of sustaining 
the velocity Vc- 

(3) Plastic failure pattern and state of stress 
that satisfy the kinematic , stress, and 
geometric boundary conditions. 

(4) Roller sinkage z measured normal to the 
original surface. 

(5) Specific rolling energy required per unit 
surface normal load and unit distance of 
travel. 

(6) Safety factor with respect to roller im- 
mobilization, which occurs when the slip 
factor S]<. = 0, for a given slope a. 

This study does not apply to the initial stages of 
rolling motion, and to conditions following roller 
immobilization. In the former case, inertial 
forces predominate and in the latter case, con- 
tinued sinkage takes place. Both cases are asso- 
ciated with unsteady conditions. This study con- 
siders the problem of soil-roller interaction only 
under uniform operating velocity conditions up to 
and including immobilization. 

B. Soil- Roller Velocity and Boundary Conditions 

As mentioned, Wong and Reece (Ref. 11) and 
Boucherie (Ref. 10) have revealed the general 
nature of the soil failure pattern associated with 
rigid driven rollers and wheels moving on level 
soft surfaces at various slip and loading conditions. 
This photographic evidence indicates (Fig. 1): 

(1) Typical leading and trailing plastic re- 
gions are formed, each, respectively, 
moving fore and aft of the advancing 
roller. 

(2) Both flow regions tend to meet at a com- 
mon point M on the roller rim surface. 

(3) Soil sinkage increases with increasing 
loads and slip. 

(4) The roller leaves no rut due to local sink- 
age. This means, after the passage of 
the roller, the soil surface fully recovers 
due to backward soil transport produced 
by the advancing roller. 

(5) The relative sizes and pattern of the soil 
leading and trailing plastic regions de- 
pend on the kinematics of soil-roller in- 
teraction. At 100% slip, the leading plas- 
tic region disappears, leaving only the 
trailing plastic zone (Fig. 3). 
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(1Z) 


In the following, the theory of sliplines is 
applied to soils for rolling with plastic flow. The 
leading and trailing plastic regions are identified 
by M(ML)EFL and M(MN)ABN, respectively (Fig. 
4). The corresponding lines M(ML)EF and M(MN) 
AB are considered to be first and second charac- 
teristic velocity discontinuity lines, respectively, 
separating the lower stationary rigid region from 
the plastic deforming one. These two character- 
istic lines meet the soil- roller interface point M 
at an angle (tt / 2 ) - <t>. It is assumed that the soil 
point M is at rest relative to the roller and that it 
pertains to the rigid stationary domain. 

The coordinates of a point i on the roller rim 
are given by (Fig. 5) 


x. = R cos [a + £.) 


(5) 


P M = tan 


- I f U z , M 


= 0 ' 

u x, M/ M 


The negative sign in Eq. (11) is consistent 
with the sign convention that a positive velocity 
component along the first slipline, when rotated 
through an angle (tt/ Z) + <(>, produces a positive 
velocity component along the second slipline (Ref. 
7) (Fig. 6b). It is of interest to note that Wong 
(Ref. 16) postulated that the trajectory of the trail- 
ing soil particle at the bifurcation point M coin- 
cides with the direction 0M of the trailing slipline 
at M. Here, the adoption of a velocity component 
Vjvf> instead of a resultant velocity as referred to 
by Wong, is due to the fact that the latter must be 
oriented at an angle <t> to the discontinuity line to 
comply with the requirements of the theory of 
plastic potential. Consequently, the trailing soil 
particle velocity resultant at M is 


z. - R sin (a* + £. ) 
1 

(6) 

T V’ 

V T = — 

M cos <p 

(13) 

The corresponding absolute velocity components 


parallel to the coordinate axes x, z are, 
respectively. 


Continuity conditions at the soil- roller 

inter- 

U . = 0)(z. - z) 

X, 1 1 

(7) 

face dictate that the radial velocity component of 
the leading soil particle at M must be equal to the 
radial velocity of the roller rim point M. Then, 
the leading soil particle radial velocity component 





(Fig. 5) is 


u . = u)(x - x. ) 

Z, 1 l' 

(8) 

V.. -r, = V. cos (a +■ §_ . - p_ .) = V w cos 
M, R M M r M M 

A 

The resultant velocity is 


M 




(14) 

/ ? ? 


where 


V. = \ u . + uf 

(9) 



1 Y X, 1 Z,1 

The angular orientation of Vj is 


(15) 



and the corresponding velocity resultant is 


Pi ''" (»;,) 

(10) 

V M ‘ V M,K ■“ [®M - <« * V * +] 

(16) 


Next, the rim boundary velocities will be re- 
lated to the soil flow in terms of the slipline veloc- 
ity components along the soil-roller interface. 

When the soil stress- strain law is derived, 
applying the concept of plastic potential, Drucker 
and Prager (Ref. 13) show that the relative par- 
ticle velocity along a velocity discontinuity line is 
oriented at an angle $ to the line. This concept 
will be used here to construct the plastic 
configuration. 


The velocity component along the first slipline at 
M is 


V M- V M C °'* < 17 > 


The angular orientation of the first slipline M(ML) 
at M is 


At the bifurcation point M, the corresponding 
trailing soil particle is subjected to a velocity 
component Vjy[ tangent at M to the discontinuity 
slipline M(MN) and equal in magnitude and direc- 
tion to the roller rim velocity V From Eqs. (9) 
and (10), for rim point i = M (Fig. 5), 


e u = e- +y - 

M M2 


(18) 


Equations (11) to (18) represent the soil velocity 
boundary conditions at the plastic bifurcation point 


M. 


V 


M 


« 


M + U z, M 


1/2 


= -V! 


M 


(ID 


The validity of the theory of plastic potential 
to define the soil strain rates and the assumption 
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regarding the boundary velocity condition (Eqs. 11 
and 12) have to be properly justified on an experi- 
mental basis. The former was adopted because of 
its rather simple application, and the latter ap- 
pears to represent a reasonable fact. 

G. Velocity Fields 


then, from Eq. (23), = 0, indicating that there 

is no leading failure zone, which is verified from 
tests (Fig. 3) and is schematically shown in Fig. 7. 
These velocity criteria produce, in general, a tan- 
gential soil velocity "jump" at M, which may be 
obtained by projecting the velocities, F,qs.(19) and 
(22), along the roller tangential direction at M. 


The plastic failure pattern configuration is 
further postulated as follows (Fig. 4): The trail- 
ing region is divided into three plastic zones, 
M(MN)N. N(MN)A, and NAB, identified as the 
active, transition, and passive zones, respectively. 
The same applies to the leading plastic region for 
zones M(ML)L, L(Ml,)E, and LEF, respectively. 

1. Active zones . As mentioned previously, 
the sliplines M(MN)AB and M(ML)EF separate the 
rigid stationary boundary from the plastic deform- 
ing one. Since a boundary at rest must have a zero 
local normal velocity component, the velocity must 
be inclined at an angle 4 to the discontinuity line. 
Along M(ML)EF, the first slipline discontinuity, 

V 1 = 0, and along M(MN)AB, the second slipline 
discontinuity, V* = 0. Since both these discon- 
tinuity lines start at rim point M, the resultant 
soil particle velocity at point M (trailing zone) is 
(Fig. 5) 


V a 

„T V M M 

V_. = r = r = “h, 

M cos 4> cos tp M 


Next, we discuss the nature of the remaining 
soil-roller interface velocity conditions and the 
plastic domain. The governing velocity equations 
that refer to the first and second characteristic 
lines, as determined by Shield (Ref, 7), are 


d V * - (V* tan 4 4 V ! sec 4) d0 


dV ' + (V 1 tan 4 + V* sec 4 ) d9' 


(24) 


(25) 


Applying Eq. (25) to the trailing zone discontinuity 
line, with V-- - 0, yields 


dV 1 + V 1 tan 4 dD' 


(19) 



where 


M 


M cos 


a M = [ (x M ' K)l + (Z M - Z)2 1 


1/2 


( 20 ) 


( 21 ) 


V “ V M exp " e ' * tan ^ ( 26 > 


which, with Eq. (19), reduces to 

V’ = 10 cos 4 r exp [(6j^ - 8' ) tan 4] (27) 


From Eq. (20), when 4 = 0, rjyj = aj^j and Ij coin- 
cides with 1 and the trailing spiral slipline trans- 
forms into circles centered at I. 


Equations (26) and (27) indicate that the velocities 
along the discontinuity line M(MN) vary exponen- 
tially and that they relate to a logarithmic spiral 
function 


The resultant velocity of the leading soil par- 
ticle at M, based on Eqs. (14) and (15), is 


r M exp t tan $ f6 vr ‘ 


M 


(28) 


V., = V „ , cos A,, = toa, , cot A = or 
MM M M M M 


where, as shown in Fig. 5, 


( 22 ) 


The coordinate positions of the trailing spiral 
pole I[ (xj , Z'j) are (Fig. 5) 


x, - x. - r cos 8 
1 M M M 


(29) 


r . . - a, cot A 
MM M 


(23) 


£1 = z, , - r. . sin 0. , 
1 M M M 


(30) 


Expressions (19) and (22) indicate the transforma- 
tion of the rim velocity V ^ into the corresponding 
soil particle velocity at the roller interface point 
M. Particularly when 


with (xj^, z |^) given by Eqs. (5) and (6) for i - M 
and 6jvf given by Eq. (18). Also, for the leading 
zone along the spiral velocity discontinuity line, 
M(ML) is from Eq. (24), 


- 0, 


M 


= M 


M 


V* - V*j exp [(6 - 0 M > tan 4 ] (31) 
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From Kqs. (17) and (ZZ). Eq. (31) reduces to 

V* - io cos 4> exp [(6 - 6^) tan 4>] (3Z) 

which applies along 

~ exp [ie - e M > tan 4>] (3 3) 


Point Mi at 4 m 1 3tt/4 remains fixed 

under all conditions of rolling without 
slip (s^ - ' ■ 0>' 

(3) There is a conformal geometric relation- 
ship of 4 m and the TqTj line's. When 
point M shifts along the roller periphery, 
the segment TqTj moves parallel to it- 
self. Thus, for increasing (or decreas- 
ing) 4 m . the segment TqTj moves closer 
to (or further away from) the roller cen- 
ter C. 


The coordinate positions of the leading spiral pole 
1^(5^, 7-i) are 


M 


' r M cos °M 


(34) 


'M 


r sin 0 ' . 
M M 


(35) 


Equations (36) to (35) yield the geometry of the 
velocity distribution along the trailing and leading 
discontinuity lines. They indicate also that the 
velocities along these discontinuities correspond 
to rotations about the poles I] and I>, respectively. 

I he geometrical character of the poles I] and 
12, in terms of the relative variations of sj^ and 
4 m- may bo graphically described with reference 
to Fig. 8 as follows: 


The locus of poles Ij for s^ - constant 
describes a series of circles with radius 
- R tan 4 1 (Fig. 9). The center of 
these circles are located along a line 
CqC ] q oriented at an angle (tt/Z) - (j> to 
the x- axi s . Along this line, the center 
marked Cp q corresponds to the circle 
with s^ - 1.0 and the center Cp - C to 

the circle for which - 0. The center 
of a generic circle corresponding to any 
0 < s^ <1.0 will be located proportion- 
ally between points CqC) q as shown in 
Fig. 9. 

(4) Given a fixed slip-factor s^. an increase 
(or reduction) of 4\4 correspondingly 
produces an increase (or reduction) of 
both spiral radial vectors ryj and ry. 

(5) When 4\1 is fixed, an increase of s]^ pro- 
duces a reduction of ry and ry. 


(11 Given a point M on the roller rim with 

4m anc * 0 - s k - 1.0, the locus _of 

all corresponding trailing spiral poles 
(Eqs. 39 and 30) defines a straight seg- 
ment TqT i oriented at an anglq (tt/Z) - tf> 
to the x-axis. The extreme points Tq and 
T ] of this segment correspond to s-^ - 0 
and s-g 1, respectively. Any interme- 
diate point between Tq and T] pertains to 
an sg such as 0 < sj^ < 1. For 4 m > tt/Z 
the poles I j and I ^ are always located in- 
side the rol_ler periphery. The quadrant 
position of Ij is given by 



A similar expression applies to 1^ for (3^ 

(Fig. 5). In general, it may be verified 

that, for s, c R sin 4«*> 
k > 

P T § 2 (37) 


(Z) On the same basis as in (1) above. Fig. 8 
indicates that the locus of all correspond- 
ing leading spiral poles 1^ (Eqs. 33 and 
34) defines a curve MM] where M and M] 
relate to Sc - 0 anti = 1. 0, respec- 

tively. In general, points M and M ] are 
always located on the roller rim. Any 
intermediate point between M and M] 
r or res pond s to such as 0 < s^ < 1. 0. 


The geometrical implication of the above 
statements relates to the basic fact that the limit- 
ing soil stress and spiral orientations (character- 
istics) are controlled by the relative position of 
poles I] and 1^. as defined by the parametric set 
4 iM and s k- I n Section V, it will be shown that 
once the equilibrium equations and stress boundary 
conditions are satisfied, the final magnitude of the 
soil stress will also relate directly to 4\( and sp 
in terms of r^i (Eq. Z0) and ry (Eq. Z3). 

Next, consider the active plastic domains be- 
tween the soil-roller interface and the discontinuity 
lines M(MN) and M(ML). which are fully described 
by a set of radial and logarithmic spiral character- 
istics with poles I[ and 1^, respectively. First, it 
will be shown that the soil-roller velocity boundary 
conditions are also satisfied. 

As mentioned at the beginning of this section, 
V* " 0 along the trailing velocity discontinuity 
slipline spiral M(MN). The corresponding inter- 
secting first sliplines are radial lines through the 
spiral pole Ij; based on Eq. (Z5), V* - 0 through- 
out the active trailing zone. This implies that at 
a point i on the soil- roller rim interface (Fig. 10), 
the resultant soil particle velocity V?" must be nor- 
mal to the radial sliplines. Then, the velocity 
component along the spiral second slipline through 
point i is 


V ' 


V. 

1 

cos (3 . 

) 


COS 4> 


(38) 


and 


V? 

i 


0 


( 3 1 1 ) 
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Pj -- P'j - e: + <t> 


Pi 


- °i + 


(40) 


with Vj obtained from Eq. (9) and p j from Eq. (101 
and 


e. 

i 



(41) 


The velocity of any soil particle (ij) located at 
a point i alone a generic spiral slipline passing 
through rim point i within the trailing plastic do- 
main (Fig. 10) is 


V.'. - V! exp f(6. - 6. .) tan <f>] (4 2) 

i i ' hi J 


and 


pertaining to the active leading and trailing zones 
have been completely defined. It has also been 
established that the velocity of each point on the 
roller rim can be transfer mod into an equivalent 
admissible velocity along the corresponding slip- 
lines that intersect the point. In particular, the 
velocity orientation of any point on the soil- roller 
rim interface, with respect to its center of in- 
stantaneous rotation I. may also he defined in 
terms of the spiral pole positions I \ (xi , ~z\ ) and 
1 2 ("*2' ~ > ) for points along the arcs MN and ML., 
respectively. 

The locations of poles 1^ and 1 2 are functions 
of the velocity boundary conditions, described by 
sk and £jyp Instead the location of rim points M 
and I. can only be obtained in connection with the 
solution of the stress equilibrium equations. 

These equations allow the determination of the ex- 
tent and configuration of the plastic domain by 
specifying a final set of slipline directions consist- 
ing of (Fig. 4) 


YF u 0 
‘.I 


(43) 


| 9 I.O’ e i.L’ 6 M’ 0 NN' 6 No} 


Following the same continuity criteria set as 
fur the leading soil particle at the bifurcation point 
M (Eq. 14), the radial velocity component of the 
soil particle at any point i along the leading soil- 
roller interface (ML). Fig. 10, is 

V. D = V. cos (a I P. - p. ) (44) 

1 , R i r i 

whe re 

V i R C0S * 

v r cos [e. - (o 44.) if] 


Using Fq. (44) and 0; + <j> - 8' + ir/2, Eq. (45) re- 
duces to 


which will be calculated in Section VI- A. 

T’he velocity field determined thus far satisfies 
the soil- roller velocity boundary conditions and the 
velocity field equations (Eqs. 24 and 25). Their 
detailed numerical evaluation is not a prerequisite 
for the solution of the limiting stress equations, 
but their existence is important to the correct 
statement of a "complete solution 1 within the con- 
text of the theory of plasticity (Ref. 12). Further- 
more, the existence of an admissible velocity field 
must also satisfy the postulate of positive rate of 
dilation (Appendix). 

Points N and I, are singular points on which 
the velocities are. multiv alued. At these points, 
the soil- roller rim interface separates itself from 
the trailing and leading traction- free soil surfaces. 


also, 


V? -- V. 
1 1 


cos (a 4 F. - p , } 

— . , ; — T TC CC)S O 

sin (o 4 |. - © ! ) 


v! 


0 


(46) 


(47) 


For a soil particle (ij) located at point j along a 
spiral slipline passing through rim point i, 

Vfj " ° X P || an ^ " OjjJ (48) 


and 


V). - 0 (49) 

ij 


The active zone velocity expressions (42), 
(43), (48), and (49) extend up to and including the 
points along the radial lines N(MN) and I, (ML). 
Thus far. the characteristic lines anil velocities 


Rim points (NN) and (LL) correspond to slip- 
lines N(MN) and L(ML), respectively, and as such 
both pertain to the soil- roller interface. Points 
(NO] and (LO) correspond to the characteristic 
lines NA and LE separating the transition from 
the passive zones and belong to the traction- free 
trailing and leading soil surfaces, respectively. 

The geometry of the characteristic lines 
M(MN) and L(Mi.) are determined by 


and 


whe re 
X MN 


X MN ' X N 


N 


cos 0 


(50) 


NN 


x, - - x . 
ML L 


L cob 6) 


LO 


N 


+ r y e xp h. °NN I tan 4>j cos 6 


(51) 

NN 

(521 
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X ML 


and 


and 


+ r 


M 


exp 


[ (C LL ' 6 M> tan *\ 



(53) 

(54) 

(55) 


with coordinates x, z obtained from Fqs. (5) and 
(6) for i = N and i - L. 


Z. Transition and passive zones . The char- 
acteristics of the transition and passive zones 
(Fig. 4) are calculated after the soil-roller equili- 
brium conditions corresponding to the active zones 
are satisfied. This will be referred to again in 
Section V-B-3 and in Part II. 
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V. SOIL- ROLLER EQUILIBRIUM 


A. Limiting Soil Stress Field 

An admissible velocity characteristic field 
has been defined that satisfies both the velocity 
boundary conditions and the governing differential 
equations for the velocities. When the velocities 
are derived from the yield stress condition accord- 
ing to the concept of plastic potential, the stress 
characteristics coincide with the characteristics 
of the velocity (Ref. 7). The problem now is to 
associate an admissible stress field that satisfies 
the equilibrium conditions along the velocity 
characteristics . 

The limiting state of stress at a point occurs 
when the Mohr circle of stress becomes tangent to 
the soil strength envelope line defined by the 
Coulomb formula t = c + <r cot <(>, where c is the 
cohesion and <f> the soil angle of internal friction. 

In terms of the point stresses T x , <r ,, and T xz> 
considering compression stress as positive, the 
Coulomb yield criterion for soils is (Ref. 17) 

(Fig. 6a) 




1/2 


- -V (tr + <T ) sin 4> 
2 x z 


c cos 4> 


(56) 


Neglecting inertia forces, the stress equilib- 
rium equations in Cartesian coordinates x, z 
satisfy the condition 


9cr 

x 

9x 


9t 


3 z. 


(57) 


9 <r 

z 

9 7 


+ 



y 


(58) 


Equations (56 - 58) represent a statically deter- 
mined problem and constitute a hyperbolic system 
having two families of characteristics as sliplines. 
Any point in the physical plane is crossed by a set 
of two sUplines, each inclined at an angle p _ it / 4 
- 4>/2 to the directions of the major principal 
stress o'} (Fig. 6). The first slipline is conven- 
tionally identified by a counterclockwise rotation 
p from the O'} principal stress direction. The 
angular direction of the first and second sliplines 
are defined by 0 and 0’ measured positive in a 
counterclockwise sense from the x and c axis in 
the physical and stress planes, respectively. The 
yield condition (Eq. 56) is satisfied if the state of 
stress is expressed by 

= pQ + sin 4> sin (29 + 4>)J - <: cot 4> (59) 

O' - p^l - sin $ sin (29 + 4> )]] - c cot <J> (60) 

r -psin<t>cos(29+((>) (61) 

xz. 


where p is the reduced mean stress parameter 
defined by 


P 


a l - 

2 sin <j> 



+ a^) + c cot 4> 2 O 


(62) 


To completely define a limiting state of stress at 
a point (x, z). it is sufficient to know the stress 
parameters p and 0. Substituting Eqs. (59 - 61) 
in Eqs. (57) and (58) and adopting s and s', first 
and second sliplines, respectively, as a new set 
of curvilinear coordinates, the partial differential 
equations of limiting equilibrium are obtained 
along the characteristics (Refs, 6 and 19): 


9P 

■57 


s? 


+ 2p tan <j> 


- 2p tan <)> 


— - v sin (0 I- <f>) (63) 

o s 

aa 

g— ; y sin (0 - 4) (64) 


Although during plastic deformation the soil den- 
sity decreases (Appendix), it is assumed in Eqs. 
(63) and (64) that y remains constant. 

B. Stresses Along Sliplines 

Using the Mohr limiting stress circle, the 
normal stress cr and shear stress r at any point 
along a slipline may be expressed by 


a = p[l + sin $ sin (20 - 2fi + d>)J - c cot (65) 
T ~ p sin 4> cos (29 - Zfl + ^) (66) 


where U is the direction of the slipline normal 
through the point under consideration. In particu 
lar, along a first slipline 


n - o' - 


(67 ) 


and along a second slipline 


n 9 + 4> 


tbS) 


Substituting Eq. (67) or (68) in Eqs. (65) and (bo), 
the general state of stress along both families ot 
characteristics are defined by 

it - p cos^ 4 - c cot <(> (t>9) 


t ±p sin 4> cos <f> (7 0) 


The plus sign in Eq . (70) indicates that the stress 
resultant at the point is rotated counterclockwise 

I 1 
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with reference to the local normal. In what fol- 
lows. the corresponding p values in Eqs. (69) and 
(70) will be defined along the sliplines L(ML), 
(ML)M. M(MN), and (MN)N, and the stress along 
these sliplines will be used to evaluate the corre- 
sponding .soil reaction forces. 

1 . Stress parameter p along straight slip- 
lines . Alternate equilibrium equations are ob- 
tained setting ds •- dx/cos G and ds' - dx/cos 8’ 
in Eqs. (6 3) and (64): 

p + 2p tan 6 d0 \ dx (tan 0 tan <j> ) (71) 

p - 2p tan 4> d6 - y dx (tan G‘ - tan <6 ) (72) 

Expressions (7 1 ) and (72) will be used to deter- 
mine the stresses along the straight sliplines 
N(.MN) and L(ML), respectively. Since and 

arfJ constant. dG - d0‘ 0. and Eqs. (71) 

and (72) reduce to 


P \1N " P NN + Y(Un G NN + tan W(x MN " X N } 


- n + A n (79 ) 

f nn 1 mn ' ” 

whe re 

Ap MN V(tan e NN + tan < t ,)(x MN " V ” ° (80) 


Both points L and N constitute singular stress 
points around which the stress are multivalued, 
each representing geometri rally a limiting first 
and second characteristic lines, respectively, 
having zero radius of curvature. When points L 
and N belong to the soil- rim interface sid(', they 
are identified as ( LL) arid (NN). If thrust' points 
belong to the t ra cti on- f re e soil surface, they are 
identified as points (LO) and (NO). To completely 
determine Eqs. (77) and (79). it is required to 
define the value of Pj^-p and Applying Eq. 

(71) to point L, 



dp y dx (tan 6 + tan 4>) 

(73) 

dp - 

- <ip tan $ dG 



dp r. y dx (tan 0 r - tan <)>) 

(74) 

1 n p - 

-Z tan <t>6 1 


Integration 

of Eqs. (73) and (74) gives, respectively, 

For point (LO), with 
obtain 

0 6 go and p - 

PLO' wo 


p ■ yltan 0 »- tan 6)x t C 

(75) 

c b : ]n 

P L o'“ lan *°LO 



p - y(tan 0' - tan (|>)x + C' 

(76) 

and 



To determine the constant of integration 
point (LL), x = X] - 9' 6j j . and p - 

at rim 
p LL , then 

In P = 
P I.O 

Z tan 4>(6 LO - 6) 


c 

P LL ' Y<tan 0 NN " tan 6,x 

L 

For rim point (LL). 

with G G j. , p 

|>j j . and 


For rim point (NN). x - Q Qjqjq, and 

P ‘ P NN * then 


P NN 4 Y ' ,an °NN 


tan d> )x 


For point (ML) along L(ML), x = and 


P ML 


P I 1 + V( tan - ta n <(il(x ML - x L ) 


P EE + Ap ML 


(77) 


whe re 

Ap ML Y ' tan e LL ' tan W,X ML - X L> 2 ° (78) 


Similarly, for trailing point (MN) along the slip- 
line N(\1N), x x vh\j and 


P LL ~ p LO cxp [ (6 LO - e EE> tan *] 18 l) 

Similarly for rim point (NN), 

P NN p NO RXp [ (f) NN - 0 NO’ tan +] f8Z) 


The. stress parameters q and PNO are deter- 
mined, considering a passive failure condition. 
Since the soil surface is free of stresses, trg - 
t q = 0, then from the Mohr circle of stress, 


e cot 4> 

P EO " p NO " 1 - “s in 4> ‘ p o 


(81) 


According to Eqs. (77) and (79), the stress param- 
eter p varies linearly with depth; consequently, 
we can operate with average stresses along 
L(ML) and N(MN) sliplines. The average stress 
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sliplines, multiply Eq . ( o 3 ) by exp [0 tan t(>] and 
Eq. (64) by exp [-6 tan 4>] and express p in inte- 
gral. form along the first and second character- 
istic lines, respectively (Ref. 6): 


p - \ exp (-1 tan <t>6 ) 


: J ' exp j( 2 tan <!>G) J 


(6 + <(>) 


ds + C (90) 


p = y t-xp 


' exp j( 2 tan <t>6)j 
X J " exp j(- 2 tan <|>6 )j - 


ds 1 + C’ (91) 


where ds and ds 1 are the elemental arc lengths 
along a first and second slipline, respectively. 
To determine the state of stress along the spiral 
sliplines, adopt as positive directions of the 
characteristic spirals the ones which relate to 
decreasing values of 6 and 0' (Fig. 11). Ex- 
pressing ds and ds 1 in terms of its radius of 
curvatures p and p’ , for the trailing M(MN] 
spiral, we get 


ds ■ - p d6 


r M exp [ <0 M “ 6) tan +] 


d6 (92) 


and, for the leading (ML)M spiral, 


ds 1 - -p' dG 1 


r M RXp [ (0 M ' 6) tan4> ] 


Replacing Eq. (92) in Eq. (90) yields the stress 
parameter pL along the leading spiral as a func- 
tion of 0: 


J * ?xp ^(30' tan <j>) cos 0' dB 1 + C 


£ | exp j( 30 ' tan 4> )j( 3 tan <j> cos 0 1 + sin 0 f ) j + C (94) 


gration is 


3 tan 4> cos + sin O' ^ 


p ML ' C 1 


+ 1 ) cos 
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Replacing Eq. (95) in Eq. (94) results in 


p F - 4 C F - C s L exp [(O' - 6^j) tan (3 tan <t cos 8' + sin 8') 


To 


obtain p^. setting in Eq. (97) S' - 0^ gives 


with 


and 


Equation (98) reduces to 


Also, with Eq. (77), 


where 


P M = P ML + C f - C s !3 tan + COS e M ■ sin e M> 


C F - 3 tan <f> cos 61, + sin 6' 
M r M M 


F, = C F - C L cJ" 
3 1 s M 


P M -- P ML + F 3 


P M - P LL + Ap ML + F 3 " P LL + F 4 


F 4 = Ap ML + F 3 

T 

Similarly, replacing Eq, (93) in Eq. (9 1 )» the stress parameter p y along the trailing spiral is 


T ^ r M 


s TT 

COS <t> 


exp [(20 + 0j^) tan 4>J J " exp ^-36 tan <j>] cos 8 d6 + C 


Integrating , 


xp [(20 H 8 M ) tan | exp [-36 tan 4> )J( - 3 


T "V r w 

p - 5 ; , — exp I ( 20 H 8 ) tan 4> ] { exp |( - 36 tan 4> )|( - 3 tan 4> cos 8 + sin 6 ) } + C 

s (9 tail <)> + 1) cos 3 <f> 1 M J 1 1 J ' 


When 8 = 6 » p g ■= p^> and the constant of integration is 


c p M - C J ( ' 3 tan + cos 6 m + sin V = Pm ' c j 


whe re 


and 


V r 


M 


(9 tan 3 <f> t 1 ) cos 3 


T T 

C 1 = C s (-3tan<(>cos 0^ + sin 6 ^ ) 


Replacing Eq. (106) in Eq. (105) results in 


pj Pj^l - cj + C F exp (6^j - 0) tan (-3 tan 4 1 cos 8 + sin 8) 


(97) 


(98) 


( 99 ) 


( 100 ) 


( 101 ) 


( 102 ) 


(103) 


(104) 


(105) 


(106) 


(107) 


( 108 ) 


(109) 
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To obtain p MN , set 6 = in Eq. (108): 


P MN " P M ' C 1 + C J eXp [ (6 M ' W tan ( ' 3 tan * cos 6 mm + Sln e ww> = + A P> 


NN NN' “ h M p MN 


( 110 ) 


where 


Ap 


MN 1 a 




c7 + cj exp [<e M - e NN ) tan <j>J (-3 tan 4> cos » NN + sin 0^^) 


( 111 ) 


With Eq. (102) and based on Eq. (100), 


MN 


p LL + F 4 + Ap MN " PLL f F 11 2 p 0 


( 112 ) 


where is given by Eq. (103), and 


F 


1 1 


F 4 + Ap MN 


(113) 


The last condition in Expression (112) insures 
that point (MN) is at or below the ground surface. 

The equations derived in this section will be 
utilized to determine the soil stresses and re- 
actions that are considered in Sections V-C and 
E, and VI- A. 

3. Transition and passive zones (Fig. 4). 

The stress characteristic lines corresponding to 
both the transition and passive zones are deter- 
mined by solving numerically the finite difference 
form of Eqs. (63) and (64), subject to their corre- 
sponding boundary stress conditions. 

The boundary conditions of the transition zones 
relate to: (1) the stress along the characteristic 
lines L(ML) and N(MN), respectively, and (2) the 
state of stress around the singular points L and N 
as defined in Section V-B-l. These stresses are 
known from the solution of the soil- roller active 
zones equilibrium equations in connection with the 
stress compatibility conditions, described in Sec- 
tion Vl-A. For the passive zones, the solution is 
obtained by extending the transition stress charac- 
teristics into the passive zones, considering the 
traction- free leading and trailing soil surfaces LF 
and NB. This solution also yields the final de- 
formed and statically correct configuration of the 
soil surfaces LF and NB (Fig. 4). 


For any solution to a given soil-roller problem, it 
must be verified that the values corresponding to 
Expression (114) are within acceptable bounds. 
Determination of these bounds is necessary not 
only to validate the solution itself, but also to sys- 
tematically initiate and objectively search for 
solutions utilizing only admissible 6 values. Ob- 
viously, the final set of Expression (114) is ob- 
tained only after satisfying equilibrium and bound- 
ary conditions. In general, the absolute hounds to 
Expression (114) result from considering the en- 
velope of values originating from (1) the maximum 
obliquity of the stress resultant as related to the 
radial direction at the soil-roller interface. (2) 
the maximum free surface slope as dictated by 
either the local roller rim surface tangent or the 
soil natural angle of repose, and (3) the stress 
compatibility requirements along the transition 
slipline zone. Each of these cases are considered 
next. 


1 . Soil- roller interface — maximum stress 
obliquity. At any point along the soil- roller inter- 
face, the shear stress t and the reduced stress 
cr + c cot <(> define a stress resultant 





(cr + c cot 



(115) 


The solution procedure and corresponding ap- 
plications will be described in detail in Part II of 
this study. 

C. Limiting Slipline Directions 

There exist definite limitations regarding the 
values that the 0 parameters may acquire along 
the soil- roller interface. It will be shown that 
the solution of the soil-roller problem entails, to 
a large extent, calculation of an admissible set of 
characteristic directions at points L, M, and N 
defined by 


The obliquity angle 6 of q r is measured positive 
for a clockwise rotation relative to the rim sur- 
face normal (radial direction) and is defined by 

6 - tan 1 ( — — — — — - 4> (• 

\<r + ccot<p/ 

When 6 - , Eq. (116) indicates that one of the 

sliplines becomes tangent to the roller rim sur- 
face at the point under consideration. In this 
case, for the trailing point (NN). with 


FLO’ 0 LL’ 6 M’ e NN’ 6 NO( 


(114) 


$ n < i 


(117) 
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- a 


e 


N'N 


+ e N ' * 


(118) 


For the leading rim point (LL), with 


fll9) 

6 LL = « + l L + (120) 


The corresponding siipline directions for points 
(N'N) and (LL) are defined by 


then each of the Eqs. (125a, b) yields two roots 
e l* e 2. which relate to two distinct rim points 
£p, respectively, as derived either from Eq, 

(118) or (120) as required. At these two rim 
points the sliplines are tangent to the roller sur- 
face. The decision on which root to adopt is 
based on the condition that at no point along the 
rim surface will the spirals cross into the rigid 
roller body. Consequently the problem is to ver- 
ify if between these two rim- slipline tangency 
points, the sliplines diverge inward to or outward 
from the roller center. The condition precluding 
the spiral sliplines to cross into the rigid roller 
body depends on the relative dimensions of the 
spiral radius of curvature to the roller radius R. 
To this effect the radius of curvature at point N 
must be 


R sin fa + 

tan ®NN ~ R cos fa I- £ n ) - Xj 


( 121 ) 


or 


P NN 


r M exp [ (6 M ~ W tan *] 3 

COS (f) 


(127) 


NN 


+ a - sin 


-R {y 'l COS 0 NN - X 1 Sin W] 

( 122 ) 


and 


R sin (a + I ) - z ? 

tan q L = ~~ 

R cos (o- f 5. ) - x ^ 


(123) 


or 


With roots 6^ 1 and 6jq 2 . such that considering 
6 N, 1 _> fl N, 2-’ it can be shown that if pjq ] < R, 
then 6 jsj 1 is the limiting direction. Thus, - 

© fq j. But if p | > R, there is no geometrical 
limitation to the position of point N along the rim 
since all sliplines diverge outwards from the 
roller rim. Similar criteria are exercised to 
determine the limiting position of rim point L 
when the spiral pole position satisfies Eq. (120). 
Here, in order to prevent the spiral from cross- 
ing into the roller, the spiral radius of curvature 
must be 


4 T + 


. - i’i OI 

Sin — (z^ cos p 


M 


exp 


[<e 


LL ‘ 0 M ) 


tan 


<t> 


LL 


2R (128) 


Adopting Eq. (118) in connection with Eq. (121) 
and Eq. (120) with Eq. (123), the limiting 0^N 
and 0j | values for 6 - 4 1 are obtained from the 
following quadratic equations: 


( x i + A) 


COS 0 NN + 2Rz l sin * cos ®NN 


sin*" if - x 2 j 0 (125a) 


-2 n, 


(~2 + 4 ) 

t (r. Z sin Z <f> - x Z j . 0 (125b) 

When the discriminant of Eqs. (125a., b) 


cos f s ’ n ^ COS ®I I 


- 2 4 . - 2 

X 1 + "l 


_2 , _2 

x 2 + * 2 , 


R sin 4> - Rj, 


(126) 


When = D, from Eq. (21), ajyj = R and, from 
Eq. (20), 


r 


M 


R 

cos <jj 


In this case, Condition (127) is always satisfied in 
the sense that the trailing spiral never diverges 
into the roller body. 

Graphically, Conditions (118) and (120) are 
defined when the pole coordinates of Ij and 
(Fig. 13) are located at or outside a circle with 
center C and radius R sin <f>, herein called the 
"4> circle. " Tangent lines are drawn from poles 
1 1 and I? which intersects the <jj circle at points T 1 , 
T 2 for pole Ii and at points Tj, T 2 for 12- Inter- 
section of these lines with the roller rim defines 
the leading points Lj, and trailing points Nj, 

N? whose angles satisfy, respectively. 

Conditions (118) and (120). When Ij or I? coin- 
cides with the if) circle, 1 = 4 l 2 ant ^ 4 n1 = 

^N2> then Eqs. (127) and (128) apply. 


When the discriminant of Eqs. (125a, b) is 
negative, then the spirals will not cross into the 
roller body. This corresponds to pole Ij or 
being at or inside the "<j> circle. " 
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An admissible value of 0j ^ must also satisfy 
the conditions of Eq. (78): 


A PML - ^ X ML - X l’ <tan 6 Ll - tan + ) 2 0 

(129) 

When < tt / 2 , (x^^ - ) > 0; therefore, 


min ®LL 2 ^ 


(130) 


When ey ■-= 0, 4^ - and X ML - X L : th - en there 
is no leading plastic zone (Figs. 3 and 7). In this 
case, Eq. (78) reduces to ApjvlL, = 0. and the 
spiral M(MN) is tangent to the roller rim at M: 


“lo ■ “'ll ' 9 L, - " li u -7 


(131) 


Also, Eq. (81) defines p ~ p - p„ as given 
by Eq. (83). ELMO 

A similar analysis with reference to Eq. (83) 
yields, for c 2 0, 


including vertical slopes as long as the critical 
height is not exceeded (Ref. 17, p. 152). 

The transition plastic zone separating the 
active from the passive zones may, in the limit, 
disappear and allow the latter two zones to merge 
side by side; therefore, 


0 LO 2 6 LL 


®NO " 6 NN 


(138) 

(139) 


In terms of slipline parameters in general, with 
p = tt/ 4 - 4>/ 2, Eqs. (134) and (135) reduce to 


a LO ' 0 LO + ^ 


y — 0 

NO NO 


- P 


(140) 

(141) 


Based on Eq. (141), Conditions (134) through (139) 
reduce, for cohesive soils (c > 0), to 


max 9 nn £ w - 


(132) 


tt + l L + !-^ 6 Lo 2e LE 


(142) 


In general, £ jvj s tt/ 2, which, from Eq. (121), 
corresponds to 


«+£ . i + u < 6 

a N 2 ^ NO' °NN 


(143) 


1/ COS ^ " z ^ 
NN t ' an (-R sin or - x, 


and, for cohesionless soils (c = 0), with Condi- 
(133) tions (136) and (137), 


2. Free soil surface — singular points (LO) 
and (NO) . The limits' of and 0jsjQ are dic- 

tated by the condition that the free surface slopes 
at points I, and N defined by and oj^jq ca n at 

most be tangent to the roller rim surface. Thus, 


tt + 4> - p 2 ®LO 2 11 ' (P + 4> ) (144) 


p - ^ s e NO £ 4> + p 


(145) 


LO 


NO 


" + * L + i 

(134) 

" + Vi 

(135) 


Also when dealing with cohesionless soils (c = 0), 
the free surface maximum slope cannot exceed 
the soil natural angle of repose. Therefore, 


for 0L O , 0’jj 

r,i ®NN' and ®r 

y[Q. Recapitulating 

Subsections ( 

3- 1 and C - 2, tl 

ie limiting slipline 


directions at points I, and M are as follows: 
For c 2 0, from Expressions (118) and (130), 


4> < e^ L < a + l L + <(> 





obtain 

Lol 

< 4> 

(136) 





" + ^N 

No| 

< 4> 

(137) 



(146) 

Based on Expressions (120), (132), and (133), we 

(147) 


U - <j> 


NN 


tan 


^ / R cos a - T j \ 
\ - R sin a - x j / 


(148) 


Obviously, if c > 0, the angle of repose has no 
significance since the soil is stable up to and 


Between Conditions (147) and (148), the lower of 
the two upper limits shown is selected. 
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Expressions (146) to (148) apply to both cohesive 
and purely frictional soils. 

Admissible values of 9 ^q anc ^ ®NO a ^ so de- 
pend on the soil characteristics. For c > 0, 


internal forces is to offset the position of the ver- 
tical soil reaction to accommodate the rolling 
torque M - instead, when P > 0 in Fq. (153), 

the leading and trailing forces make up for the dif- 
ference in allowing for P to become equilibrated. 


and 


a 4 


■ + 2 
’L + 2 


9 LO £ 0 LL 


(149) 


Taking moments with respect to roller axle 
center C, with positive torques measured counter- 
clockwise, moment equilibrium requires 


a + 


NO 


NN 


(150) 


✓M = 


M. L 4 

k 


M. 


M r = 


(155) 


And for c = 0, it can be proved that Conditions 
(149) and (150) also apply but subject to the follow- 
ing limitations: 

ti + ^ - p > 0^ o 2 ir - (p + <f>) (151) 


where Mj,. is the applied axle torque, and mI‘ and 
correspond to leading (I.) and trailing (T) mo- 
ments produced by the mobilized soil strength 
and corresponding soil weight. The above forces 
and moments relate exclusively to the active zones 
L(ML)M and M(MN)N. Expressions (153), (154), 
and (155) are further developed as follows: 


and 


p - 4> s e NQ < 4 . + p (152) 


For lunar locomotion (c > 0), then Expressions 
(149) and (150) will generally apply. 

Expressions (146) to (152) are incorporated 
in the computer program to operate within a com- 
patible range of stress parameters d and 0. 

D. Equilibrium Equations 


With the positive direction of forces the same 
as the positive x, z coordinate directions, and 
with applied vertical axle loads W and pull force P 
parallel to the terrain slope, horizontal equilib- 
rium conditions require (Fig. 11) 


4 4 P cos a - 0 (153) 

and, for vertical equilibrium, 


Lz = W, L 4 wj + W + P sin cv — 0 (154) 

k K 


where 11^ , , Wg are horizontal and ver- 

tical forces corresponding to leading (L) and 
trailing (T ) plastic active zones adjoining the soil- 
roller interface arc I.N. Sub-index k (=1,2, • • • ) 
identifies the possible existence of more than one 
soil- roller solution for a fixed slip value and 
various 

A driven roller moving over horizontal or 
sloping terrains operates under self- propulsion 
conditions when it transports only its axle weight 
W (pull force P “ 0). Under self- propulsion 
there is no net soil thrust development nor soil 
resistance to motion since the leading and trailing 
forces balance each other (IlJ" 1 I IlJ' - 0). The 
net effect of the interplay of these self- equilibrated 


2x - 11^ + h7 4 + H T 

k , s k , s R , p k , p 


4 P co s a - 0 


(1561 


+ W. L 4 V, r , L + W, T 4 w T 
k,s k.p k,Y k , s k,p 


+ W. + W + P sin a - 0 

k, Y 


(157) 


Lm = M, L + M, L 4 M, L + M. T 4 Ml T 

k,s k,p k,Y k f s k,p 


4 M, T - M. = 0 

k, y k 


(158) 


The sub- index s corresponds to forces (or 
moments) due to stresses along the spiral slip- 
line M(ML) and M(MN). The sub- index p indi- 
cates forces (or moments) due to stresses along 
the straight slipline I, (MI.) and M(MN) as derived 
from both the transition and passive zones. 
and y are soil weights within the confines or 
the leading and trailing active zones M(ML)L and 
M{MN)N, respectively. Mj^y and MjJ'y are 
moments due to W-j^y and w’jk.y j respectively. 

To evaluate these forces and moments it is neces- 
sary to determine the nature and extent of both 
active plastic domains (Fig. 11), their asso- 
ciated stresses and corresponding soil weight 
participating with the roller motion. 


E. Soil Reactions 


The horizontal and vertical force components 
due to stresses acting along slipline L[ML) are 
(Figs. I l and 12) 
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H. 


i - — ■ qi — ni l LI L L i 

(tr sin 0 ' - t cos 0 ' . ) — ■ n i — - b 


k.p “ ' L LL L LL cos 0' 


LL 


(159) 


where 


X MN = X 1 + r M exp [< 0 M - 0 NN ) *] COS 9 


NN 

(169) 


W F - - ( 0 - cos 9' + t sin 0 ' ) —^^- 57 — — b 

k, p L LL L LL cos 


x N = R cos (o- + § N ) 


(170) 


(160) 


where b is the roller width, and 

‘ xp [ |0 LL ‘ V tan *] cos 6 


, , = x, + r, , ei 

ML 2 M 


LL 

(161) 


Substituting Eqs. (169) and (170) in Eqs. (167) and 
(168) and ordering terms results in 


H k,p = C 15 p LL sin <0 NN + + F 13( 0 NN ) (17U 


x - R cos (o- + | ) 
1—1 


(162) k,p 


Replacing Eqs. (161) and (162) in Eqs. (159) and 
(160) and ordering terms results in 


w v „ ■= - C 15 [ p N COS (0 NN + ♦> - C ° S ®Nn] 

(172) 


where 


H k. P = C 14 p LL Sin ( 9 Ll - + ] + F 12 ( 0 LL ) (163) 


C.- = ~b (x. ... - x ) (173) 

15 cos 0 MN N 


W 


k.p 


- c m[pl cos ( 0 Ll - - set? cos 0 , ll] 


F 13 (e NN } " C 1 5 I ( F 1 1 -^T^) sin (0 NN +4,) 


(164) 


whe re 


sin 8 , 


sin <(> NN 


(174) 


' 1 4 : b cos°^ U ML- X J (165) 


F 12 < 0 LL ) = C 14 


AP 


ML . , , 

2 sin <0- LL - 9) 


Expr e ssions (163), (164), (17 1), and (172) will be 
incorporated in equilibrium Eqs, (156), (157), and 
(158), The horizontal and vertical force compo- 
nents on the elemental arc ds along the spiral 
<ML)M are (Fig. 11) 


— - — T- Sin O' _ 
sin 9 LL 


(166) 


T T 

Derivation of the force components Hj^p, on 

the trailing first slipline N(MN) is basically similar 
to the procedure adopted for the leading slipline 
L(ML). The horizontal and vertical force compo- 
nents from stress acting along N(MN) are (Fig. 12) 


„T X MN ' X N , 

H k,p = sm0 NN + " COS 0 NN cose" b 


dH, 


" - -b[cr L cos (0' - 9) - t l sin ( 6 1 - 9) ds 

: , s L s s j 


(175) 


dW F s - -b sin ( 0 1 - 9 ) + cos ( 0 1 - <t> )j ds 


(176) 


(167) 


,..T „ _ . „ X MN ' X N , 

W k, P = °"t COS 9 NN ' T sm 0 NN ) cos TT ~ b 


NN 


( 168 ) 


where cr|^ and correspond to a and t as 
given by Eqs. (69) and (70) respectively, with 
p - pL (Eq. 97). 

Similarly, the horizontal and vertical compo- 
nents on the elemental arc ds 1 along the spiral 
M(MN) are 
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dll/ bfuTsin 0' - tJ cos O') ds' 

k . s * ^ 


dW 


.1 


k, s 


-bttr^ cos 9 ! +tJ sin 0 1 ) d s : 


U 77 | 


<178) 


where trj and t ^ correspond to (r and r as 
given by Eqs. 169) and (70), respectively, with 
p pj (Eq. 108). 

The total force components are derived by 
replacing the cr.r stresses in Eqs. 1175) to (178) 
and integrating within the corresponding spiral 
limits. Thus, 


H 


k . s M 


W 


i b r . . 

k , s M 


'’ 9 m 


LL 


0 M 


EL 


K<- 


CDS (0' - 26) 


cos (9’ - 6) 


i exp (9' 


r L 

[ p s S 


sin (0' - 26) 


sin 6 


sin ( 0 ' - 6) 


I exp j^< 0 ' 


9 , , ) tan 6 d9' 

M 


0j^j ) t an 6j d0 ' 


(179) 


(180) 


NN , 


H. 


k, s ' br M 


W, b.r, , 

k, s M 


p cos e 


TiTT* cos (e + 6) J exp [ (0 m ’ 0) tan 6 J de 


(181) 


M 


NN r . 


jp^ sin 0 - ^ . [ e sin (9 + 6(J exp ^(Ojyf - 9) tan bj dG 


(182) 


M 


Replacing Eq. (97) in Eqs. (179) and (280) and Eq. (108) in Eqs. (181) and (182), integrating, and 
arranging terms, the soil reactions s and Wv g along the sliplines are obtained, i or the leading 

zone, the solutions are 


and 


<s - C s[pme' ^L'V^i^LL'V 


(183) 


w r. s C 4- p ML* F 7 <e LE’ e M )+F 8 <9 LE’ 0 M 


(184) 


wh f 


C B - br M cos 6 


(185) 


- c - jsin O' - exp [(O' - 0 ' ) tan 6] sin 0' 

sin <J» [ M 1 [ LL M Li LL 1 

+ c] J jsin (9^ + 6) - exp j^0' LL - 0^) tan sin (0' LL +6)j 


tan * 4 Mr + sin (29 m + *>] + 1 tan + 4 exp [ 2(0 ll - 9 M ) tan 6 ] 


x [nsr 1 sin I20 'it + 4>) ] f | cos (20 m + 4,1 - ex p [ 2 ( 0 Ll - 0 M ) tan *] 


X cos (20j j + 6)j 


(186) 
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F 2 - Sin (0 M f ' eX P [ t0 LL * 0 M ) tan +] Sin ,0 LL + 4,1 
F-. -- cos (9^ + 4>) - exp - ej^) tan 4>j cos (6j L + <(>) 


8 sin 4> 


! COS 9 M - (>xp ( le LL - e M> tan +] c ' os e LLi - C 'r | cos ,e M + ‘ f ’ ) 


- exp ^(0'j^ - 0J^) tan 4>j cos <0' IjL f <t>) J + ^ tan <j> C F jeos (29^ t 4>) 


- exp [ 2(0 'll - e M ) tan +] cos , 20 Lf + *>| - ■ sin (20 m + 4>) 


-f ex P [ 2(0 'lL - V tan * 


J x ■ sin i + 4 > )1 

sin L.L J 


For the trailing zone, 


T 

" C 7 

k, s T 


[ p M F 5 (9 NN-V +F 6 (e NN’V 


and 


W k,s " C 7 [ P M 1 9 ( °NN’ 0 M^ + 1 10 (9 NN' 0 m'] 

where 

C„ - br„ , cos cb 
i M 

F 5 ‘ exp [ (0 M - 9 NN ) tan ♦] sin (0 NN - + > ' Sin (0 M - +> 

F 6 t - sT^TJ | e> l-’ [ (0 M - *W tan +] Sin 0 NN - Sm %! 

- Cj jexp [(0 M - 0 NN ) tan *] sin (9^ - ) - sin <0 M - 4>)j 

- ? c s exp [ 2(0 m - 0 nn ) tan +] | 3 tan 41 [- TiTT - * + sin (2S nn - +) ] + c ° s <20 nn - 

+ i c I | 3 tan * [- + sin (Z9 m - +>] + cos (29 m - *>i 

r 9 = exp [ (0 M ' e NN> tan *] cos (0 NN ' 40 ' 003 (0 M ’ 40 

F 10 ; ilF? | eXp [ (6 M ■ 9 NN> tan *] COS 0 NN - C ° S 0 Ml 

+ cj jexp [<9 m - © NN ) tan 4>] cos (9 NN - 4> ) - cos (9 M - 4>)j 

+ i C J ex P [ 2 < 9 m ' °nn ) tan ♦J [ 3 tan * cos (20 nn - 4,1 ' 7777^ + sin (29 nn - ?) ] 

■K[s tan 6 cos (20 M - 4») - ^Xj + sin <26 M - * )] 


(187) 

(188) 


(189) 

(190) 

(19D 

( 192 ) 

(193) 


(194) 

(195) 


(196) 
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Expressions (183) to (196) will be incorpo- 
rated in Eqs. (156), (157), and (158) for the de- 
tailed study of equilibrium in Section VI- A. 


For the trailing zone M(MN)N, the soil 
weight is (Fig. 1 1 ) 


F. Body Forces 

Soil body forces W^y are obtained by inte- 
grating the soil volume contained within the lead- 
ing and trailing active zone L(ML)M and M(MN)N 
as shown in Fig. 11. 


W k T ^ = "yb j^area IjM(MN) 


area circular sector I j MnJ (205) 


For the leading zone, 


r.v = Yb [ : 


area ML)M 


Area IjMflVlN) 


]r 2 de 


- area circular sector I,LM (197) 


where r is given by Eq. (28). 


Area - 


ir 2 d9' (198) 


Area IjMN - (area triangle IjMN + area 
circular segment with arc MN) = or 


= i r r sin (0-. - 0 . .) + A (207) 
y 2 M e M NN c 


where r is given by Eq. (33). 

Area I 2 LM = (area triangle I^LM + area 
circular segment with arc LM) = Cy or 


where 


A C = 7 r %^n (20s) 


C Y - i r U 7 c sin(0 LL- V + A c (199) 


! 1 N - R cos (a I- - 0 NN ) 


A c 2 R Jr ' ' Sm ( T " W] 


r c , 1 Z L - -R cos {c + S L - e' LL ) 


- a T COS - P T ) 


n - hA 

P T = COS (id 


+ a L cos (p L - 


a T = J x 1 + Z 1 


(3^ - tan 


-1/ “2’ 


4 N = 6 NN * a+s ' n ~ l [i lz l cOS Vn - X 1 Sin ®NN ) . 


ljt, _2 
a L ^ J X 2 + Z 2 


Integrating Eq. (206), and with Eq. (207), Eq. 
(205 ) becomes 


Integrating Eq. (198), and with Eq. (199), Eq. (197) 
becomes 


( * 

jrd£* | t,xp [ 2(0 ll - 9 m ) tan 4 


with CL given by Eq. (199). 


I 4 tan <)> 


Jexp [ 2 ( 0 M - 0 NN ) tan *] -i] - C^j 


where is given by Eq. (207). 

Equations (204) and (213) wi U be used in Eqs. 
(15h), (157), and (158) for the detailed study of 
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soil-roller equilibrium in Section VI- A and also 
to determine the driving torque M in Section V-G. 

G. Moments 


( 218 ) 
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For the trailing zone M(MN)N, 


M 


M k,> = 


,fl NN 


4 r 2 jx^ + -j r cos e| d8 - 


(224) 


where r is given by Eq. (28) . and 


M ' t- r_ .r sin<6.. -0- I . T )(x.,+x. 1 4x.)+Ax 
c 6 M c M N N M M 1 c c 


A c is given by Eq. (208). and 


. 3 I^M ' 


sin 


COS V 


1 


4 \ 2 / — - \ 2 

N ' " 3 6m - - Sin ( ^M - i 


N' 


Integrating Eq. (224), 


M, - Vb 
k. y 


- 2 _ 
r M X 
4 tan 


+ !-P [ 2( °M - °NN ) tan +] - 1 j 


M 


' M 


3(9 tan 4> + 1) 


jexp ^(» M - e NN ) tan 4>] ( 3 tan * cos 6^ - sinG NN ) + (sin 6 m . 3 tan <j> cos 0 M )J 


(225) 


The driving torque will be used in Section 
VI- B to determine the roller driving power 
requirements. 

H. Soil-Roller Interface Stresses 

Once the pattern and extent of the active plastic 
domain zone are derived, the soil-roller interface 
stress at a generic rim point i are (Eqs. 65 
and 66) 


<r. 

i 

- p. |I + sin <f> sin 4 hj - c cot 43 

(226) 


t. r p. sin <b cos 4 * . 
1 *1 ' 1 

(227) 

wh ere 


- ze - 2 ft + 4> 

i 1 1 

(228) 

with 


n i “ + h 

(229) 


For the leading zone, 

®i * + h ’ sin ' '[i( S 2 COS 6 i ' *2 Si ” 6 i)] 


and for the trailing zone, 

6 . = a + 4 . - sin (z, cos 0 . - x. sin 0 .) 
i i LEI i 1 i . 

which results in 


F -1 - - s in’ hi ( z , cos 0 ! - Jf , s in 6 ! )1 t ir - <() 

l [R 2 l 2 l J 

(230) 

and 

+ 4 - - sin’ ^i(zj cos 6, - sin 6,)j + tf> (231) 


The parameter p along the slipline M(ML) is 
derived frorrvEq. (9?), and for the slipline M(MN) 
from Eq. (108). On the leading zone, along a 
radial slipline which intersects the roller rim at 
a point i(4j. > 4m) ar >d also the sliplino M(ML) at 
point (Mj) is 


pJ^PMi+AP^ <232) 

where ! J Mj corresponds to Eq. (97) for S' - 8.', 
and 

Apf" = \(x i - x Mi )(tan 6| + tan <f>) (233) 
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with 


Mi 


<2 + r M exp 


fa 


- & M> tan * 


cos e! (234) 


and xj from Eq. (5) and from Eq. (23). 

Similar criteria apply to the trailing zone 
( £, j < tt/2). Along the radial slipline intersecting 
the roller rim at point i and the slipline M(MN) 
at point (iM) , 


T 

Pi 


P: 


iM + *Pi 


(235) 


where pT^ corresponds to Eq. (108) for 6 = 0j, 
and 

T 

^P i = y(x i - x.j^Xtan G. - tan 6) (236) 

with 

x iM r X 1 + r M eX P K * 6 i J tan +] COS B i {23?) 

and xj from Eq. (5) and ry from Eq. (20). 

The orientation 6j of the soil- roller interface 
stress resultant q r is given by Eq. (122) in con- 
nection with Eqs. (226) and (227). 
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VI. SOLUTION OF BASIC EQUATIONS 


A. Basic Equations 


and 


It was shown in Section IV - C that the geometry 
and extent of the active and transition zones re- 
late exclusively to the slipline directions at points 
L, M, and N, as defined by Expression (114) 

(Fig. 4). Also, if the set (Expression (114) is 
known, a velocity slipline pattern can be unambigu- 
ously built that satisfies the roller velocity con- 
ditions. To initiate a solution of a given soil- 
roller problem, as postulated in Section 1V-A, a 
set of parameters and is selected first. 
These parameters define unique values of 0^ 

(Eq. 12) and 0{v[ (Eq. 18) and also determine the 
position of poles I] (Eqs. 29 and 30) and I 2 (Eqs. 

34 and 35). 

To evaluate the four remaining unknown 
parameters of Expression (114), there are avail- 
able four basic simultaneous equations. Two of 
them originate from satisfying the horizontal and 
vertical equilibrium conditions, as given by Eqs. 
(156) and (157). For completeness, Eqs. (156) 
and (157) are repeated as Eqs. (238) and (239): 


C 2 = C 8 ' F 2 + C 1 4 Sln (0 LL - + C 7 • K 5 


+ Cj 5 sin (0 NN + 4>1 


Since, in Eq. (240a), ^ > pg ? 0, Gj and G 2 

must have the same sign. 

From Eqs. (81) and (83), 


P LL = P 0 ° XP | f0 LO ' 8 Ll’ * an *1 (240b) 


After equating (240a) and (240b), 


G, 


0 ' - 8 ' + 


In 


LO LL 2 t an 4> p q G ^ 


(241) 


EX = H, L + H T + H, L + H, X + P cos a 
k, s k, s k, p k, p 


F H f9 LO’ G EL' 9 M’ °Nn' ° 


Also, it is known from Eqs. (82) and (83) that 
P 0 exp [(0 N>I - 0 NQ ) tan <i>j (242) 


NN 


(238) According to Eq. (79). 


EZ = W. L + W, L +■ 
k , s k, p 


W 




T T 

W, + W. 
k, s k, p 


P NN P MN " Ap MN 


(243) 


+ W 


T 

k.V 


+ W I P sin a 


F w ( 0 Lo’ 


°Ll- 0 M> W = ° 


Alter equating (242) and (243), 


0 NO 6 NN 


1 


2 tan 


In 


P MN ' Ap MN 


(244) 


(239) 


Substituting in Eq. (238) the corresponding 11^ 
force components given by Eqs. (163). (171), (183), 
and (190), after a short trans formation, results in 


I.L 


Gj ~ F L (0 i,L' G NN* 


(240a) 


where 


G 1 " ' C 8 <F 2 ’ Ap Ml, + F l } ‘ F 1 2 


' C 7' F 4 ' F 5 + F 6>- F 1 3 " PcOS ° 


Thus, we arrive at a system of four equations 
(Eqs. 239, 240a, 241, and 244) in the unknown 
parameters 0 lO> ®LL' ®NN. 0NO- These equa- 
tions are solved by iteration using a computer 
program. Admissible values of 0V r are first 
selected within the corresponding bounds stated in 
Section V-C. Subsequently, ' s determined 

from Eq. (239). ^ Substituting GlL an d ®NN ' n Eq. 
(241), a unique 6 value is determined which sat- 
isfies simultaneously the horizontal and vertical 
equilibrium equations, With ©EO. an d ®NN 

known. ©NO i s determined from Eq. (244). During 
the solution process, the computer prugram verifies 
the fact that the determined angular parameters 0 
comprise only admissible values; otherwise, a 
new 6 L L is selected or a new case (sj<, £m) 
initiated, as required. Equations (238) and (239) 
express the necessary conditions under which an 
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Equilibrium solution is admissible with regard to 
the soil-roller active plastic domain. It is of 
significance to point out that the existence of the 
set of values, Expression (114), which satisfies 
the equilibrium and the stress compatibility equa- 
tions, does not necessarily imply that a solution 
to the problem has been found, unless the follow- 
ing conditions are also simultaneously met: 

(1) The rate of dilation, as mentioned in 
Section IV-C must he positive throughout 
the plastic domain. (This is verified in 
the Appendix. ) 

(2) At no point outside the plastic regions 
shall the stress exceed yield. 

(3) The free-suvface points F and B (Fig. 4), 
which belong to both the rigid and plastic 
domains, must also be aligned with the 
original surface slope a\ thus, 


Z F ~ Z B 

tan a - (245) 

X F " X B 

which assumes that, after the roller 
passage, the original trailing surface 
slope is not significantly altered. This 
was verified experimentally by Wong and 
Reece (Ref. II) for roller tests on level 
surfaces. Here, this condition is assumed 
to prevail also for slopes. 


is proven when, in addition, the solution of the 
two active zones can be appropriately coupled to 
the neighboring transition and passive fields, 
which satisfy Eqs. (63) and (64) and the remaining 
boundary conditions. Here, we advance the fact 
that the active plastic field solution can be ex- 
tended up to and including the traction -free sur- 
faces LF and NB (Fig. 4), if the input data is con- 
sistent. In Part II it will be shown that the 
completion of the plastic field also yields the 
statically correct deformed free surface and that 
Condition (245) can, in general, be satisfied, 

B . Sp e c i f i c Energy Dissipation 

The roller moves parallel to the original soil 
surface with uniform velocity V c - loRpS^ where 
R e is the effective rolling radius. For a rigid 
cylindrical roller R e = R. The axle traverses a 
distance H per unit time, 

L - ojR s, (246) 

6 K. 

The total soil thrust parallel to the original sur- 
face is 

T P + W sin a (247) 

with P the drawbar pull force and W = the ap- 
plied axle load. The load component normal to the 
original surface is 


Although conditions (1) and (2) do not bear 
directly on the solution process of the equations, 
they are fundamental to describing the require- 
ments for the completeness of a solution. Con- 
dition (2) is not specifically verified, but may be 
considered satisfied if no sharp corners of rigid 
material develop within the rigid plastic boundary. 
In fact, it may be proved that the soil state of 
stress at point M does not exceed yield and, in 
general, it may be assumed that the rigid material- 
can support the plastic deforming body. The 
method of extending the plastic stress field into 
the rigid domain, as studied by Shield (Ref. 20) 
and by Cox, et al. (Ref. 21), can also be applied 
to this problem. 


Regarding Eq. (245), the coordinates of F 
and B are obtained, satisfying the boundary con- 
ditions on the leading and trailing traction-free 
surface slope. The positions of points F and B are 
satisfactorily approximated based on the following 
analysis. It was noted in Sections V-D and 
V-E that the horizontal and vertical stresses along 
L(MIj) (Fig. 11) determine the reaction forces 
and , respectively. Similarly, the 

stresses along M(MN) define the force components 
IIj^ and These leading and trailing forces 

have a special significance with regard to the 
validity of any solution of Eqs. (238) and (239). If 
a solution is found, it has to be verified that the 
mentioned forces p and p can be sustained 

at the corresponding transition and passive plastic 
domains. In other words, it must be verified that 
all plastic zones satisfy simultaneously the basic 
requirements of local and overall equilibrium. In 
this context, a solution to Eqs. (238) and (239) 
represents only a necessary condition. Sufficiency 


N = W cos a 


(248) 


The total energy input Ej^ per unit time due to an 
applied torque M at the roller axle is consumed by 
the thrust force energy output F.y and the soil- 
roller energy dissipation Eg/j^: 


E M ' E T + L S/R 


(240) 


where 


E. . = Mu 
M 


(2501 


and 


E t = TL 


(25 11 


In general, the soil energy dissipation per unit 
travel length Eq. (246) and per unit normal load 
Eq. (248) will be defined as the soil- roller speci- 
fic energy dissipation coefficient*. 


J S/R 

ML 


WR s 


_M ( 

3 k cos * ' V 


W cos a 


+ tan ol 


(252.) 

Equation (252) accounts for the soil and tire 
deformation and soi 1-wheel interface friction energy 
losses due to slip. It represents a general energy ex- 
pression and define s the pe rformance of any rolling, 
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power-driven device; the equation applies to both 
rigid and flexible rollers or wheels. If the roller 
is flexible. E will express the specific energy- 
dissipation produced by both the soil and the roll- 
ing device. Care must be taken to properly mea- 
sure or calculate the effective turning radius R e . 
For the rigid roller, the values of M and S]^ in 
Eq . (252) are determined front the problem solu- 
tion (Section VI-A). 

When a 0, Eq. (252) reduces to 


F _.-M . 

J 0 " WR s, W 
e k 


(253) 


Equation (253) was used by Leflaive (Ref. 22) to 
analyze test results of driven rigid and flexible 
wheels on horizontal soil surfaces. The specific 
energy ( Eq . 252) shows two terms. One is the 
specific torque energy input at the axle per unit 
norm a) load and unit travel distance. 


E 


M 


M 


WR s, cos o 
e k 


(254) 


The other is the specific thrust energy output per 
unit normal load and unit travel distance. 


E 


T 


P 

W cos a 


I tan cr 


(255) 


Only when a =; 0 docs the specific thrust energy 
equal the pull/load ratio Eq- = P/W; otherwise, 
for a > 0, Eq. (255) refers to the specific thrust 
energy. Equations (254) and (255) can be used to 
evaluate the performance of power -driven vehicles , 
providing the relative wheel slip factors and axle 
load distribution are known. The wheel thrust 
efficiency is given in general by 


n% 


_P 

W 


+ sin a 


M 

WR 

c 


s, X 100 
k 


(256) 


In practice, the specific power consumption per 
kilometer of travel along a straight line on a slope 
a 2 0 is given by 

<**» 

where Em is defined by Eq. (254) for > 0 and 
W is the applied axle load in Newtons. 

C. Rigid Roller Sinkage 

The roller sinkage z is measured perpendi- 
cular to the- original surface. Once the leading 
and trailing points F and 13 are determined, satis- 
fying the basic equilibrium equations, it is veri- 
fied if points F and B are aligned on a slope a 
(Fig. 14). To this effect, the coordinates of F and 
B yield 


tan a ' 


where d ' is the direction BF for a trial solution. 
Only when a' a does the solution found corre- 
spond to the given problem. Then, with 


t = 


+ Xg tan a 


and 


n - l cos a 


the sinkage is 


z 


R - n 


(257) 


The derivation of Eq. (257) presumes the 
condition stated in connection with Eq. (245). 

D. Mobility Safety Factors 

It was shown in Section Vl-A that the soil- 
roller mobility problem can be solved satisfying 
the velocity and limiting equilibrium equations that 
are subjected to the corresponding boundary con- 
ditions. The solution determines the operational 
slip factor and torque M for a given axle load- 
ing P and W. It was also shown (Section I) that, 
when - 0, the roller becomes immobilized 
(Vq - 0). From a mobility safety standpoint, 
given the soil conditions (c, <f>, V, a) and a fixed 
set of operational loads, pull Pq and weight Wq, 
it is required to determine the corresponding 
maximum load W max (or P max ) which, in combina- 
tion with Pq and Wq, produces immobilization of 
the roller. Thus, there exist two basic loading 
conditions capable of immobilizing the roller: 

(1) increasing only the pull force from Pq to 
Pmax. and (2) increasing only the axle weight 
from Wq to W-j^g^. Consequently, two types of 
mobility safety factors (SF) can be defined, de- 
pending on the ultimate cause that stops the roller. 

The first definition of SF is 


SF 


P + W_ sin a 
max 0 

Pq + Wq sin a 


(258) 


where the numerator indicates that the maximum 
soil thrust Tjy,^ is reached by incrementing Po, 
pull load, to Pmax- setting sfc = 0. The denom- 
inator Tq corresponds to the roller operating 
thrust for 0 < sj^ < 1. 0. In this case, when 
Pmax r Pq, SF - 1 and the roller would stop 
due to excessive pull. 

The second definition of SF corresponds to 


SF 


I 1 . + W sin a 
0 max 


Pq + Wq si n a 


T 

max 



(259) 
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where, as in Eq. (258), the T rnax corresponds to 
s k = 0 by incrementing the axle weight Wq to 
W max - When W max = Wq, SF = 1 and the roller 
stops due to excessive weight. In essence, both 
Eqs. (258) and (259) relate to the soil maximum 
thrust capacity developed for - 0. 

For the general case of sell propulsion; 

Pq = 0, the roller propels its own weight Wq. 
Under this condition: 

(1) If a > 0, Eq. (258) reduces to 


SF 


P 

max 


+ W n 


sin a 


Wq sin a 


(260) 


In Eq. (260), if for s k = 0 it is deter- 
mined that P max = 0, then SF = 1. 

This condition defines a limiting roller 
slope angle climbing capability, 
a - <r max , for self-propulsion. 

(2) If a = 0 , then T max = 0, and Eq. (259) 
has no practical significance, since on a 
level terrain self-propulsion is unrelated 
to soil thrust. This conclusion derives 
from the fact that, under the action of a 
vertical load, there is no net soil thrust 
mobilized. Under this condition, the 
leading and trailing soil reactions are 
balanced (Eq. 153) (Fig. 11): 


II 


L 

k 


+ H, 


0 


Therefore, in this case, the SF refers 
exclusively to the maximum vertical soil 
load capacity for s k = 0. Thus, from 
Eq. (259), 


W 

SF = — (261) 
0 

Equation (261) is the factor of safety 
versus immobilization valid only for self- 
propulsion on level terrains (a - 0). 

For Pq > 0 and n > 0, the applicable SF 
definition corresponds to Eqs. (258) and (259), as 
specified. 

Given = 0 and the computer program 

determines connected with an operational 

pull load Pq. On the same basis, given an opera- 
tional load Wq, it is possible to determine the 
corresponding P max , which immobilizes the 
roller. 

E, Applications 


dimensions (R,B). Two additional input param- 
eters, 5^ and slip factor s k , arc also required 
and are entered by means of the computer’s tele- 
printer, The user, with a minimum of experi- 
ence and iterating on and s k values, can 
determine a number of possible admissible solu- 
tions. In this report, results obtained satisfy 
the equilibrium Eqs. (238) and (239) (active zones) 
and the limiting conditions stated in V-C-3. It is 
expected that from all solutions which are found 
at least one will satisfy Condition (245). (Refer 
to Section VI-A. ) It is also assumed that if more 
than one solution satisfies Condition (245), then 
the one with the minimum torque energy (Eq. 254) 
would represent the actual wheel performance. 

Thus far, Condition (245) is not incorporated 
in the SWIP program. This will be done only 
after the plastic field is extended into the transi- 
tion and passive regions, utilizing a program sub- 
routine (Part II). 

The SWIP program outputs the following: 

(1) Stress parameters p, 9 and soil-wheel 
rim interface normal and tangential 
(shear) stresses (psi) at points L, M, 
and N (Fig. 4) . 

(2) Geometric parameters. Coordinates of 
the center of instantaneous rotation I, 
the Reading and trailing spiral poles Ij 
and I;>, and the radial angular directions 

and where the soil detaches from 
the wheel rim at the leading and trailing 
edges, respectively. 

(3) The total and partial, vertical and hori- 
zontal soil reaction forces considering 
the effect of soil weight on both the lead- 
ing and trailing regions (to an accuracy 
of 0. 5%), The partial moments a:id total 
required driving torque. The specific 
energy input Ej^ (Eq. 254) and specific 
thrust energy Ep (=E-p) (Eq. 255). 

The program outputs the nature of any incom- 
patibility which may arise from either the equilib- 
rium Eqs. (238) and (239) or from the limiting 
conditions stated in Section V-C-3. By this means, 
the user can, on an interactive basis, select 
appropriate new sj. values to bring about a 

solution. 

In what follows, the SWIP program is applied 
to estimate: 


(1) Driven rigid wheel performance tests car- 
ried out on horizontal terrains under con- 
trolled slip. The wheel axle is subjected 
to a vertical load Wq and a pull force Pq 
parallel to the undisturbed soil surface. 


The foregoing soil-roller analysis was pro- 
grammed in Fortran II for use with the IBM 1620 
computer. In the following, it is assumed that the 
soil- roller model developed also applies to a 
finite-width roller (wheel) as long as the predom- 
inant soil failure mode occurs in the fore -aft 
direction rather than in the lateral direction. The 
soil-wheel interaction performance (SWIP) pro- 
gram input data is; wheel axle loads (W,P), sur- 
face slope a, soil properties (Y, 4>, c )< and wheel 


(2) Driven rigid wheel slope climbing per- 
formance. The wheel axis is subjected, 
to a total vertical load W = (W^ + P$) ' 
acting on a slope angle defined by a ~ 
tan"^ (Pq/Wq), where Wq and Pq are 
loads corresponding to the horizontal test 
conditions defined in ( 1) , above. The 
purpose of this slope climbing calculation 
is to verify if there is theoretically any 
performance difference when equivalent 
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wheel normal and pull loads act either on 
a horizontal or a sloping terrain. 2 

(3) Lunar roving vehicles (LRVs) on the 

assumption of driven rigid wheels rolling 
on a level lunar surface. This is applied 
particularly to the Apollo and Lunokhod- 1 
vehicles. 

Table 1.0 provides a summary of the above- 
mentioned applications indicating typical wheel 
loads, dimensions, and soil properties to be used 
in connection with the given application. 

Test case 1, Table 1. 0, was performed by the 
Waterways Experiment Station (WES), Vicksburg, 
Mississippi, under controlled 25% slip (s^ - 0,75) 
(Ref. 18). The SWIP program applied to this test 
condition produced the results shown in Tables 1. 1 
through 1.3a, which correspond to 4 m - 99, 101, 
and 102 deg. For intermediate values of 4 m> such 
as 99 deg £ 4 m £ 102 deg, there exists an infinity 
of solutions satisfying Eqs. (238) and (239). The 
indicated results correspond only to the bounding 
values pertaining to a given set of 4 m anc * s k = 

0. 75. intermediate solutions were also obtained 
but unfortunately lack of space precludes their 
inclusion. As mentioned, these results have to 
also satisfy Condition (245). This condition will 
eliminate all those cases which do not meet the 
boundary requirements referred to in Section VI-A. 
It was also found that for sk - 0. 75 and 4 m : 
deg and 4 m - 107 deg, there are no other com- 
patible solutions, thus indicating the fact that the 
operational range on 4 m ' s bounded and if a solu- 
tion exists satisfying Condition (245), it must lie 
within the results given in Tables 1. 1 through 1. 3a. 

The measured torque was M - 600 lb-in. and 
the results indicate that this value is appropriately 
bounded by the program output as shown. Addi- 
tional rigid soil wheel tests were also checked 
using the SWIP program. Particularly, in the 
case of W 108 lb, P - 30 lb, and 25% slip 
(Ref. 18), the measured torque was M - 720 
lb-in. and also was satisfactorily approximated 
and bounded. 

In general, concerning the mobilizable soil 
strength, it is typically assumed that the same c, 
<|>soil parameters apply to both the trailing and 
leading regions. Any divergency between the lead- 
ing and trailing limiting values of c and <j> to be 
used is concerned with the question of how the soil 
parameters c and <f> are modified due to the dis- 
turbance produced by the passage of the wheel's 
leading edge or by other wheels along the same 
track. For instance, with reference to the total 
torque, a difference in <|> values does not appear to 
be very sensitive, as seen in Tables 1. I — 1. 3a 
for <(> - 4 2. 3 deg and Tables 1.4 -1.4a for <j> = 35 
deg. However, in order to satisfy Condition (245) 
and thereby arriving at a complete solution, it may 
be necessary to resort to different c, <j> values for 
the leading ami trailing zones. 

Regarding the wheel slope performance, 

Table 1.0, case 2, the computer results are 
shown in Tables 2. 1 through 2.4a. First, it was 
found that there are no compatible slope solutions 


for a 25% slip as it occurs for o = 0 test. This 
indicates that there is no analogy which relates 
equivalent loading between horizontal and sloping 
tests. Second, the slip performance for self- 
propuision (P = 0) is a minimum when a ~ 0 and 
slip increases for increasing slope angle. When 
both 4 m and are varied, as shown, solutions 
will be found between the values indicated. For 
self-propulsion conditions, the horizontal leading 
and trailing soil reactions are equal and of oppo- 
site direction. The slope climbing energy re- 
quirements are generally higher than for a - 0 
due to the combination of larger torques and wheel 
slippage. 

A hypothetical application of the SWIP pro- 
gram to the Apollo LRV flexible wheels is given 
in case 3, Table 1.0. on (he assumption of rigid 
wheel behavior. Results shown in Tables 3. 1 to 
3.2a represent self-propulsion conditions on a 
level lunar soil surface. Results indicate that the 
wheel operating range for this case is within 10".» 
to 15% slip (sk 0.90 to 0. 85). Calculations 
also indicate that the wheels could not operate at 
20% slip for a 0. The soil -wheel rim interface 
stress level is rather low (<3/4 psi), and energy 
requirements are not unlike the expected mobility 
performance for on-earth operation. 

Case 4 (Table 1. 0) represents an application 
of the SWIP program to the Lunokhod- 1 to investi- 
gate its mobility performance for a 0 . To this 
effect, lunar soil properties similar to those 
applied to the Apollo L0.V (case 3, Table 1.0) are 
considered. Under self-propulsion, it is assumed 
the wheel load is approximately W - 35 lb (15.6 
kg) (Ref. 25). The wheel radius scales roughly 
R _ 10.0 in. (25.4 cm) and width B - 6 1/2 in. 
(15.35 cm ). Th e wheel i s as sumed to operate as 
a rigid finite-width roller. A pattern of grousers 
and a metal mesh covers the wheel rim which, on 
ground contact, confines a soil layer of an approxi- 
mate thickness equivalent to the projecting grouser 
lugs. This condition insures the soil-wheel inter- 
face mechanical properties are at least equivalent 
to the lunar soil strength, thus eliminating any 
uncertainty connected with the soil-wheel inter- 
face adhesion. Since the Lunokhod- 1 mobility 
performance is independent of the wheel's surface 
material, appropriate correlations of lunar soil 
properties can be made utilizing its mobility per- 
formance records in a lunar traverse. 


The operational performance of the Lunokhod -1 
is shown in Tables 4. 1 and 4. la. Results refer to 
4 m - 112 deg for 20% slip (s^ - 0. 80). It. is noted 
that no compatible solutions were found for 4 m ' 

1 10 deg and 0. 80. 


A review of the results shown for the Apollo 
and Lunokhod- 1 vehicles indicates that after 
extending the plastic field up to the t raction -f ree 
surfaces, Condition (245) may be satisfied and a 
complete solution defined. Further applications 
of this program are planned to estimate the limit- 
ing wheel slope climbing performance conditions 
after incorporating Condition (245) (Part II). 


In this connection, vehicle mobility tests on slopes indicate that the vehicle performance degrades with 
increasing slope angles. Test results indicate that slope tests cannot strictly lie simulated by equiva- 
lent wheel axle loading performed on horizontal terrains (Ref. 1). 
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VII. CONCLUSIONS 


(1) A comprehensive theory tor the solution 
of the soil- roller interaction problem has 
been presented. This solution is appli- 
cable to power-driven rollers moving on 
horizontal or sloping soft soil surfaces 
under quite general conditions of terrain 
slope angles, soil properties (cohesion, 
friction), and loading conditions including 
gravitational effects. 

(2) In this study, Part I, the method of solu- 
tion satisfies both the roller velocity 
(slip conditions) and equilibrium require- 
ments within the active zones (Fig. 4). 
The solution was programmed for com- 
puter use. The nature of the developed 
soil-wheel interaction performance 
(SWIP) program is that it only outputs 
bounding values of wheel performance 
parameters. In Part II, it will be shown 
that these bounds can be narrowed further 
and that, from the point of view of the 
theory of plasticity, complete solutions 
can be obtained which satisfy overall 
equilibrium, velocity, and boundary con- 
ditions, which include both the transition 
and passive zones (Fig. 4). 

(3) It is considered that a finite-width roller 
also represents, on a first approxima- 
tion, the performance of a rigid wheel, 
which must be verified by tests. Limited 
application of the theory to rigid wheel 
tests on level terrains indicates that 
experimental results compare favorably 
with theoretical predictions. Experi- 
ments have to be performed considering 
mobility on level and sloping soil sur- 
faces, taking into consideration the 
underlying concepts of the theory, as 
formulated, particularly with regard to 

(Fig. 4); 

(a) Soil-wheel failure pattern on slopes. 

(b) The existence and shifting of the bi- 
furcation point M separating the 
leading from the trailing plastic zone 
along the soil-wheel interface. 

(c) The position of the leading and trail- 
ing edge detachment points I. and N, 
re spectively. 


(d) Laboratory determination and prac- 
tical use of the soil parameters c 
and cj> as related to potential soil dis- 
turbance affecting wheel performance. 

These are just a few of the many items 
which must be considered before accept- 
ing this or any other theory. 

(4) The limiting soil- roller interface radial 
and tangential (s hea r) stresses were 
defined and it was found that the obliquity 
angle of the resultant interface stresses 
with respect to the radial directions 
varies along the roller rim. 

(5) A general normalized energy Expression 
(252) was derived which is of practical 
use for evaluating and correlating wheel 
(vehicle) test results for slopes (a > 0). 
The nature of Expression (252) and 
limited application of the theory (Tables 
1.1 — 1.3a and Table s 2.1 — 2.4a) indi- 
cate that the wheel thrust, torque, and 
efficiency performances relate to the 
particular slope a. This result points to 
the fact that wheel tests using equivalent 
normal and pull forces on horizontal and 
sloping terrains do not represent similar 
loading systems since the state of stress 
and limiting equilibrium conditions of a 
soil slope and a level terrain are different. 
Thus, horizontal wheel tests based on 
equivalent loadings cannot be used to pre- 
dict wheel slope climbing performance. 

(6) A safety factor (SF) concept against wheel 
immobilization is introduced which is 
applicable to any driven rigid or flexible 
wheel for varying loads and slopes. This 
SF concept sets the framework for the 
study of mobility as a basic mechanical 
process whereby a safety number can be 
assigned to each of the commonly used 
wheel efficiency performance parameters. 

(7) Regarding the validity of the theory, par- 
ticularly its reliability, the proposed 
method of solution has to be more exten- 
sively evaluated by applying the computer 
program to a wider range of mobility con- 
ditions, slope angles, load combinations, 
soil properties, and wheel slip values. 
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VIII. RECOMMENDATIONS 


It is recommended that: 

(1) The theory and computer program devel- 
oped for driven rigid rollers (wheels) on 
soil slopes be extended to also include 
towed rigid rollers (wheels). 

(2) Both the driven and towed rigid wheel 
solutions, referred to in (1) above, be 
generalized to consider flexible driven 
and towed wheels on soil slopes, thus 
covering the whole spectrum of potential 
wheel operations as may be applicable to 
different mobility modes on planetary 
surfaces. 

(3) The solutions mentioned in (1) and (2) 
above for single wheels be coupled to 
consider the mechanical interaction 


tween the wheels of a vehicle. Since each 
wheel of a vehicle system is subjected to 
varying loads, wheel slips, torque, ter- 
rain slopes, and soil properties, predic- 
tion of vehicle performance requires 
knowledge of coupled wheel mechanical 
behavior. This program will assist in 
(a) modeling vehicle-terrain interaction; 
vehicle design configuration, safety factor 
against immobilization, and power re- 
quirements; (b) defining planetary vehicle 
operation mode s; route selection, deci- 
sion risks, and data rate requirements; 
and (c) interpreting mobility operations 
and test results. 

(4) I'he results of the theory be verified and 
validated by implementing a comprehen- 
sive soil-wheel interaction testing pro- 
gram. 
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APPENDIX 

POSITIVE RATE OF DILATION 


A. Introduction 


Drucker and Prager (Ref. 13) applied the con- 
cept of plastic potential to Eq. (56) and derived 
the stres s - strain laws connected with the rigid 
perfectly plastic material. On this basis, the 
axial plane strain rates are 


3u 

• x 

6 x 3x 


\ 

2 


j^sin 


<j> - sin (26 + <t>) 


(A- 1 ) 


Next Eq. (A-2) will be determined in connection 
with the strain rates derived from Eqs. (A-4) and 
(A-5). 

B. Trailing Zone (£ ; s ir/2 ) 


For a generic trailing point (ij) (Fig. 10), 
setting 6.. = B in Eqs. (42) and (43), 


V! . 
1 J 


V’ = 


V'. 

1 



- 9) tan 


(A-6) 


9u 

Z 

dz 


^2 1 sin <j) + sin (26 + 4>)| 


(A-2) 


v*. = V* = o 


(A-7) 


where Visa positive factor of proportionability, 
in general a function of time and position. For 
steady state X - X (x, z). The rate of dilation based 
on Eqs. (A-l) and (A-2) is 

A = 4 e = X sin <*> 2 0 (A-3) 


From Eq. (A-3), for <j> > 0, A > 0. Shield (Ref. 
7) expressed the velocity components Ux,u z of a 
point at failure in terms of the slipline velocity 
components V" ! and V as follows: 


and substituting Eqs. ( A-6 ) and (A- 7) in Eqs. (A-4) 
and (A-5) results in 


,,, sin 0 sm 6 

u = -V' r - -V. v exp 

x cos <p l cos 4> 


j^(e i - 6) tan 


(A- 8) 


u = V' - V 

z COS 


! C - ° - S ^ exp ["(6. - 6) tan 4>1 
1 cos <t> L 1 j 


(A-9) 


u 

X 


V* cos (e + <j>) - V' sin e 
cos <)> 


(A-4) 


With x, z coordinates of point (ij), 0;j 


°J 


6-. 


u 

z 


V' 1 ' sin (6 + <j>) + V 1 cos 6 
cos 4> 


(A-5) 



(A- 1 0) 


5 i 
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C. Leading Zone (£j 3- £^1 


9 x x 


exp j^( 0 . - 0 ) tan t)>j (tan 4 > sin 8 


- cos 0 )L 


(A- 11) 


For a generic leading point (ij) (Fig. 10), 
setting 0jj ~ 8 in Eqs. (48) and (49). 

Vvt - V* -- Vr exp J^(e - A ) tan <j>j 


(A-18) 


0Z Z 


V.'. - V' - 0 
i.l 


(A- 19) 


v i r 1 ae 

exp ( 0 . - 0 ) tan 4> (tan A cos 8 + sin 8 ) -r — 

cos 4> L 1 J dz 


Replacing Eqs. ( A- 1 8 ) and ( A- 1 9 ) in Eqs . (A- 4) 
and (A- 5), with 8 + 'b -■ 0 1 + tt/2, 


(A- 12) 


From Eq. (A- 10), with x - Xj = rj exp [(0j - 0) 
tan t|>] cos 0 , 


98 sin 8 
^ exp 


[(6 - 9 .) 


) tan 4> (A' I 3) 


u. :t - cos (0 4 . <K ) 

x cos cp 


V . sin 0 ' 


exp jjO - 8 j ) tan j> 


(A- 20) 


— = - ^ exp [(6 - 6 .) tan <j> (A- 14) 

1 t- 


Substituting Eqs. (A- 13) and (A- 14) in Eqs. (A- 11) 
and (A- 12), respectively, Eq. (A-2) reduces to 


u - ~t sin (0 + 4>) 

z cos ® 


VrCOS 6' _ -| 

“77: :: * eX P [ (e - V tan<t, j 

(A - 2 1 ) 


A. , - _ — — tan A (A- 15) 

i) r. cos ^ r . 


Introducing Eq . (38) in Eq. (A- 10), and with 
Eq. (40), 


With (x, z.) coordinates of (ij), 


- 1 / 7 z 2 
0 . = tan I 

1 \ x - x 2 


(A- 22) 


V . V . t an d> 

• i j ~ « 

A. . - «- tan - jzl — — jri — ; — T\ - A. 

li r. cos p. r. cos (p. - 0 . + <p) i 

J 1 ' 1 l ' r i i 

(A- 16) 


du 

X 

9x *x 


ty exp ^(0 - 6 . ) tan cf>j 

cos ij (_ 


tan 4> sin 0' - cos 0' 


V. tan 4 

A. = — , 2 0 (A- 17) 

i r. sin (0, - p^) 


9z 7 , 


Equations (A- 15), (A- lb), and (A- 17) indi- 
cate that under steady state conditions the dilation 
rate A— at a point (ij) reduces to the dilation Aj of 
a point i on the soil-roller interface. Also for 
4 > ■ 0 . Ajj - 0 , representing the incompressibility 

condition of a Tresca material with c = shear 
yield stress. 

To satisfy Eq. (A - IT), sin (0 i - p; ) must be 
greater than zero; this condition in general holds 
true as may be verified graphically or analytically 
in most practical cases. 


exp |"(e - 8.) tan 4 >j ^ 

cos ^ [_ 


tan rjj cos 6 ’ - s in G ' 


From Eqs. (A-23) and (A- 24) with x - x£ 
iq exp [(0 - 0 j ) tan <j>] cos B ' , 

88 sin 0 ' .. 


- - S — sr- — exp |"( 8 ; - 0 ) tan <(>J (A- 25) 
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exp 


(A- 26) 


36 _ cos 6 ' 
9z r j 


fe i - e) tan 


In particular, for rim point i = M, 6 - pj^j 

and Eqs. (A- 16) and (A-28) yield, with Eqs. (20) 
and (23), 


After substituting Eqs. (A-25) and (A- 26) in Eqs. 
(A-23) and (A-24), respectively, Eq. (A-2) re- 
duces to 


A T - V M tan 4> 

M r. , cos <t 
M 



tan cf> 


(A- 29) 


V5= 

A.. =. = - — -r tan (A- 27) 

Ij r cos <p ^ 



“ — cot A tan rf> 
r M 



tan 6 


(A- 30? 


and with Eq. (46), 


V. cos (a + 4 . - p.) 

A * . -■ ~ • t \ t n 1 1 tan — A . 

i j r sinfa + s. - 0 !) 1 

i 1 1 


(A- 28) 


To satisfy Eq. (A- 2 ), \a + f T - p ; - tt/2, which is 
satisfied for 0 < s]( < 1 . 0 . 

As for the trailing zone, Eqs. (A-27 ) and (A- 28) 
also indicate that the dilation rate at point (ij) 
reduces to the dilation rate Aj at point i on the 
soil-roller rim interface. 


M 



Consequently, as expected, since the state of stress 
along the soil-roller interface is uniform and 
continuous, the soil dilation rate at point M is the 
same when approaching M along (ML)M or along 
(MN)M(Fig. 4). Also, since A^ is proportional to 
(rad/s), it would be of interest to verify to what 
extent the theory of plastic potential (Kef. 13) is 
applicable to the prediction of soil deformation ui - 
der conditions of steady- state motion. It is know 
that for continued straining under unsteady condi- 
tions the dilation predictions far exceed the mea- 
sured increments (Ref. 23). 
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NOMENCLATURE 


velocity arm 
B or b roller width 

C roller center 

c cohesion 

ds,ds' elemental arc lengths along firstand 

second sliplines 

E soil-roller specific ene rgy d is sipa- 
tion coefficient 

energy input per unit time due to 
moment M at roller axis 

Es/p soil-roller energy dissipation 

Ep or Ep specific thrust force energy output 

H horizontal force soil reaction 

I center of instantaneous rotation 

Ipl^ leading and trailing spiral poles 

i point on roller rim 

ij point of intersection of i and j 

j point along slipline 

L distance parallel to slope n 

M moment 

N normal load 

P-! : roller load per unit width 

P pull force on roller axle, parallel to 

slope (P - bP !,! > 

P power consumption per kilometer 

p stress parameter 

q r stress resultant ( soil- roller 

interface) 

R roller radius 

R g effective rolling radius 


R^ radius of 1 circle" (Fig. 13) 

r,, r. leading and trailing spiral radial 

vectors of point i 

SF safety factor related to P 

’ max 

SF safety factor related to W 

7 max 

s slip % 

s, s' first and second slipline curvilinear 
coordinates 

slip factor = 1 - s 

T thrust force 

u . , u . velocity components of point i in x 
■ ’’ and z directions 

V translational velocity of roller 

center C 

V. absolute velocity of point i 

V*, V.' velocity components of point i along 
first and second sliplines 

W roller axle weight, vertical reaction 

x, z Cartesian coordinates 

z roller sinkage 

a terrain slope 

p-p, |3^ angular orientation of I j , 13 with 

~ reference to x-axis (Fig. 5) 

V soil unit volume weight 

A increment 

A rate of dilation 

angle as defined in Eq. (15) 

6 obliquity angle of q r with respect to 
a radial roller direction 

€ x’*z axial strain rates, x and z directions 

0 stress parameter 


36 
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0. , 0! 

angular orientations of first and 


thrust efficiency 

i l 

second sliplines at point i 

n 

angular orientation of a vector 

X 

positive factor of proportionability 


normal to a surface 

P 

~ it/ 4 - *(j/ 2 

10 

angular velocity, rad/s 

6i 

angular orientation of rim point i 




f Fig. 2) 

Supe r sc r ipts 



angular orientation of rim point M 
(Bifurcation point) (Fig. 4), deg 

L 

leading 

P, P 1 

spiral radii of curvature for first and 
second sliplines 

T 

trailing 



Subscripts 


P i 

angular orientation of velocity V. of 




rim point i 

k 

see Section V-D 

<x 

normal stress 

P 

pa s sive 

T 

shear stress 

s 

spiral 

4> 

soil friction angle 

T 

transition 
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Table 1.0. Soil- wheel input data for SWIP program application 




Loads 

r 

Wheel 

Soil 




Case 

Appli- 

cation 

W, lb 

P. lb 



<y , 
deg 

V. , 

lb/ in. 

<t>. 

deg 

C, 

psi 

Slip s, 
% 

Torque 

M, 

lb-in. 

R eniarks 

i 

Horizontal 

test 3 

94. 0 

35. 0 

13. 95 

12. 0 

0 

0. 0584 b 

42. 3 b 

0. 06 b 

25 

600 

Tables 1. i- 
1. 3a 



94.0 

35. 0 

13. 95 

12. 0 

0 

0. 0584 

35. 0 

0.06 

See re- 
marks 

See re- 
marks 

Tables 1. 4- 

1.4a 

2 

Slope 

100. 3 

0 

1 3. 95 

12. 0 

20. 4 

0. 0584 

42. 3 

0. 06 

See re- 
marks 

See re- 
marks 

Tables 2. I- 
2. 4a 

3 

Apollo LRV 

60. 0 

0 

16. 00 

10. 0 

0 

0. 01 

33. 0 

0. 05 

See re- 
marks 

See re- 
marks 

Tables 3. 1- 
3. 2a 

4 

Lunokhod- 1 

35. 0 

0 

10. 00 

6. 5 

0 

0. 01 

33. 0 

0. 05 

See re- 
marks 

See re- 
marks 

Tables 4. 1 
and 4. la 

a Reference 18. 
'■’Reference 24. 
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Table 1.1. Horizontal test for = 99 and = 0. 75 (upper bound) 

I m Pin a r a m P r 5 R S * * * 

>■ , l» t X I ( ; - 1 , SK = <44.0000 35.0000 90.9499 , /ftOi) 

ft L L, H'\ , G/U ■»■'»■! , t>H I , C = 0.00)0 .0 9 >‘4 43.1999 .OftOli 

K , K - ] 3, 99 00 ) 2. 1)! HI;, 


***s i rfss rara^p m*r s*** 

Till’ TA (LHP ) = 13«.RH3 o THh r A ( L-_ 9 > - 1I)U.4?1 4 i'Hf- f A ( Mp ) = Vb^Si)(i 

T H I" T ft ( •' 1 ) = Hfl.l*iMH 1'Hn r A ( rNJ.M ) = 3/.797U I I R" T A ( MU ) = 7.27/0 

PI I - . /U/9 L>l_l. = .70/9 PM'. - l.l 1 / n 

PI*' s .9936 PMN = .694 1 PNM = .40/*) 

I*-' ip MAI. S T R 5 S S F S AT L*rt f w 17 .2/14 .994 9 .2? 19 

TANOFf'iT T A I, .5IKF5.SP.5 AT L»M,N= . ? VH) .07/0 .2034 


i C P ARAMPTf-R S*«* 


CPmT'tH HP I'MSfAMT AM»J.)US RO f A T 1 1 )M 

XH = 

0. oouu 

2H = 

10.40/9 

trailing spiral phlf — 

xp L = 

-3.il 30 

/Hl = 

H .4 132 

LhAHlMf. SPIRAL PI )Lk 

X PZ = 

-3. OH/ 9 

/P/ = 

12 . / 905 

Hi PURL ATI mm point 

xw = 

1 H/Z 

2M = 

l 3 . V 7 M / 

LPAniftir; fruit 

XL = 

-3.HU40 

7L - 

13.4/11 

U-AniKi Spiral curmni-JAT KS 

X'1L = 

“4.7131 

/ Ml. = 

1 H . 1 290 

TRAILING fih;f 

X M = 

4 . 3 M / ';» 

ZN = 

l 3 . 2430 

TRAILING spiral cmiir ui mat f.S — 

XMN“ 

0.5 3/ '{ 

Zmn = 

14 .HZ 2 5 


XI (I. )= 1 OS.H/OR XI(N)= 7 1 * OH 9 / 

*** VLR T [ CAL R PACT t UN*** 


i-K Pi.= — 20 .OH// WKGL= <+.5/ LI WKSL = -39.44/5 IiiTaL LlA'ilMl. = -56.5535 

WK*>T = 1.7050 WKGT = 6.P974 WK5>T = -47.051H TOTAL T ^ A [ L I -3M.4H44 

* * =:• HH K I 7 l 1 M T A L khACFIi In* * * 


H < P l_ = / 7. OH 09 HKL = -9 .0/25 Till A i. L PAD I MO = IH.uSh'3 

HKPT= -13.1*595 MKT = -34.H4RH UiTaL TRAILING* -53.05H3 


* * * M OM F N T S * * c 

MGL = -1R.274/ MSI. = 270*7107 MP L =-345.1012 liHAl. LEAPING = -HO. 0030 
MGT = 19.0000 MST = 423.4217 MPI = 195.04/0 TOTAL THAU. INb= ft 3H. 0.304 


T 1 1 T A l. THRfjllh = 501.3007 LH-ImCHHS = 70.2/93 K (j-MHFKS 

FM,fP, = .5000 .3/23 
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Table 1. la. Horizontal test for = 99 and = 0. 75 (lower hound) 





‘•fi [ 'U>:li J A * Af- 

u T i- 3 3 ::i * 


^ » L> » * 1 ( ' 

) ♦ S 

= 

9ft .(lO'io 

3 ft , onno 

9ft . 

l\l. |M- \i\ t (\ A 

i ■■ 1 /) , 

i ; ‘H 1 , 

(. - o,ni)i)o 

• o ft ft ft 

4 3.1 999 

k t K = 




j 3 . 9ft Oil 

j / . UtJiJU 


>. f ; s- kAkaoKT F« S~- >?t 

THH A ! I.iib ) 167.1776 ThI|'a(LLH)= ]io.?-ti4 fi If T A (m P ) = 13 . 35 OH 

TH(-TMM) 1|I|. 1SIW THI- r A ( I'H'U = 2U.66K7 Tin f n I i'.ill ) = LK.nl l rt 

U|i = . 7074 PIJ. - .‘>770 HMt. - .4404 

M*-i = .‘MOO »hrt - . 'tHH'l = .7137 

nMIJHAI. STKFSSFS A ( - .7664 .07/1 .0671 

f/U'Oi-N T ! ftl SIKFSSFS A I .7416 .hilt .1171 


A I -il- IK 1C PAKAMt-THi.S «*#* 


on ih: rie i u.s i'am i ai-ihius « r i t a T 1 « — 

XD - 

0. uuou 

7H = 

1 1 ■ . > >2ft 

1 PA 1 [ IMf- 

. S P I ft A L P < 1 L ft 

X PL = 

- 3 . L L 3 h 

/.ft L = 

H . 4 i 3 2 

1 P/iiUN(, 

SPIRAL PULL 

XP2 = 

“3. 6H2 9 

/H2 = 

1 2. 7906 

rtIPilKCATTM'M 

XM - 


n - 

13 . 7 7B2 

i>ad i-'ir, 

MH'.i 

X'. = 

- 3 . 9 04 7 

7.L = 

1 5.3923 

h AD I N(; 

Si } I ft 41. CMHK!) I 'J AT 1- S 

X^L = 

-ft . h7 3'3 

= 

.1 ft . 7 33 ft 

T k A I i. I M r; 

- 1 '1 f , i 

XM = 

h, Hfih 

/ M = 

1 2. f 390 

T K A T l T 

tl UiRi) ( MA 1 ftS 

XMN = 

1 0.4 6 74 

/ ift I\1 = 

1 3 . 4 '( ft 9 


X r { l_ ) = 10k. 7647 XIliM)= 60.6146 

**»Vd-<T ! Cal K FAC 1 1 dn#» * 


W< PL - 
UK PT = 

- ] ft . h i ) ] q 
2 .^648 

i*J KG I. - 
ujKGT = 

ft ■ A ft ft 3 
l 7 • ft 7 ft 7 

WKSL 
wK S T 

= -4 / • ? 1 4 L 
= -S2.D3ftO 

mi a l 
t n i a l 

l Pad In = 
r K A 1 1. 1 to G = 

-47.349/ 
-3ft. ft 9 39 



I 7H’M f Al 1 

RfcAC 1 

I 1 IN *# - 




HK PL - 

hkpt = 

31 .937? 

-9 . 

HKL = 
hKT = 

-6. 770 1 
- 67 . 7 66H 



T 1 1 T A L 
TH1 AL 

L K A o I ML) = 

TKA 1'. 1 1611 = 

2 1 a? ft 7 ft 
-62 • 2ft / ft 




* vvi^i M i 1 NTS** 





ivif,| = 

M GT = 

-73.31 Hh 
ft7 .n^n i 

-3SL = 
MS T = 

i> 1 0 .0 Is 9 

rtKI. 
13 P T 

= -43 7. 79«j>; 
= 144.4ft9^ 

Till a L 
1 0 1 A L 

LF Ail I 'M ( t = 
I'KA 1 1. I Pile 

-2U/.?lft4 

/lft*ft293 

T OT M_ 

TMKUlUh = 

SlHJ 

.1174 1. It- 

I AlCHP 

5 = 7 0. 

6141 K(. 

'-Mt-T GPS 


ft At f ft P , 

* 


= 

. ft!7M 

. 3 /7 3 
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Table 1. Z. Horizontal test for = 101 and s, = 0. 75 (upper bound) 

<=** [ aipiit P AKAMF PiR !>»»# 

l»l,SK = 44.0000 33.01)110 10(1.4444 .7300 

AL P4A , 0 AMMa , PH I . L = n.noon . 0*104 43.1444 .0400 

0 . B - 1 3. 4000 1 2 . UOuo 


*w*s trfss j ar/h-'f r fr s «* 

THtTA(l,n°l= 130.ASS4 1 H I- I' A ( 1 .'. P ) = H/.HK74 lhFrA(MP)= -(4.4000 

THf-TOlM = 00 , ?H OH THFT* ( NO ) = 4l.0o7() f H F T A (MO ) = 14.43(7 


PH = . 20 23 

PM = 1.{H9 3 

P L l. = 

P MM = 

.8229 
.8 223 

I J M l__ - 

PlMN = 

1 . 04 96 
. 3« 1 1 


mi ip mai. sihfssfs at . i.,m,n 
1 ANHFNTTAL .STkFSSFS a 1 L, 

, H , IM= 

.1430 

-.Olio 

1 • 1 4 4 h 
. 6SlU 

• ■i 6l ) 1 
. 3 6 34 


not to i c p ao Ann r i-o Sv* * 




CFMTFR U4 \ i\jS r AMT Ai\iri.HJ S 

Rn r at n in 

XK = 

0 * UOUO 

Z* = 

10.4623 

TK A I LING 3 P Ik A I, P (1 L F 


XP I = 

- 3. 0342 

ZPl = 

7 .9629 

L F A H T M f ; SPIRAL PDL^ 


Xl^ - 

-4.414/ 

7 P2 = 

1 2 . 24 96 

H I Fl IRC AT T f iM Pui^f 


XM = 

-2.6hl / 

Zm = 

1 3 .OS 3b, 

Lf-AIHNG FDGh 


XL = 

- 4 . 3 / rt 0 

ZL = 

1 3. 24 3 1 

LFADtMr, SPIRAL CIIMRIUNATFS 

XM | = 

“4. 22 97 

Zml- 

1 7.23 73 

that iHir, foul 


XN = 

2 . 3469 

ZM = 

13.7311 

TkAlLJIxin SPIRAL COOK n 1 NAT FS — 

XMi\| = 

4 . 3 99 u 

7. F i N “ 

1 3 .9383 

X ] d J = ] nn . 29 08 

XI ( N ) = 80, 

. M 4b 





*«*VhRT I C A L RLACT I IIN**» 


uK P L = 

-23.3242 

i,i KOI. = 

3. IZ'i'i 

WKS L = 

-31 . 2H73 

total 

LfAlllMlr = 

-3 1. 4 86 2 

wkpt* 

- 1 » 6 h (3 1 

wKGT = 

5. SO 13 

WKST = 

-46.41 63 

TOT AL 

Trail Im(1= 

-42.SH12 



W^HflR I 7 DMT A L ; 

RFACT I OIM'I- *<* 




H K P L = 

20. 11)00 

HKL = 

-1 1.0933 



1 OTai. 

L FA 01 no = 

9.0017 

HKPT = 

- lb) . H 108 

H K T = 

-27 . I 90H 



TO 1 AL 

IRA 1 L INU = 

-44.001 l 




# M | H 3 N T S * 'X 4 





m.(;l = 

- 13 .9624 

MSL = 

2R4. 7330 

MPL =■ 

-207. 9697 

TOTAL 

LFAIllMG = 

39. 8229 

MGT = 

4 . H709 

MST = 

31 1.46L3 

MPT = 

24 3 . 06 Si 

1 O f AL 

T RA 1 L I Nb = 

3 6 3 « 4 0 L 4 

TOTAL 

TOR OOF = 

b2i 

.2243 L6- 

I NChFS 

= 86. 

1 04 0 <( 

.-met f kn 


FM,FP 



— 

.6333 

.3 r?3 
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Table 1. 2a. Horizontal test for - 101 and - 0.75 (lower bound) 

N p 1 1 T PARAMFT frs*s* 

w,p,xi (m),sk = OA.nnnn 34.onno ldo.vvvv .74nu 

Al. PHA ,RAWMA , PH I f C = 17.0000 .05 HA 47.1999 .0M)(l 

«,« = 13.4500 17.0000 


SS# S 1 RFSS .’MAMURS*?* 

THFTA ( l.flP )= 137.447H THl= TM U. P ) = 44.0444 IHHA(MP)- 79.4Hi>H 

THFTA(M) = H4.7H0H THFTA(NN) = 34.4334 |HFTA(MII)= 27.4447 

PI] = .2074 P|_>. = .4242 PMl_ = .94n4 

PM = , 42H3 PMIM = .4524 HNM = .2544 

M U P M A L SlRrbSLS AT = .1594 1.0156 .104 1 

TAMOFMTTAL 5T4 F S SO S Al i.,M,N= .1440 .61 < <L .1416 


sssr.ijonCTRtC parameter ssss 


CFNTLR 0- ImSI'ANI AMFfHIS 

RO TAT MIN- 

XH = 

0. oouu 

ZH = 

10.462S 

TRAILING SPIRAL PUL' 


X P 1 = 

■ 3. 0 3 42 

ZPI = 

7 . 0620 

L F AD ! MG SPIRAL POLL 


X P 2 — 

-4.4L4 { 

z p?~ 

12 • 24 06 

KU-iJRCttT TON PI II N T 


XM = 

-2. 601 V 

/m - 

1 3 .69*36 

l_ FAN INC, FOGF 


XL * 

-4 . 6 63 6 

ZL = 

13. 1H24 

L FAIT I Air. SPfKAL f.lllWDl MA TFS 

XML = 

-6. 364 3 

Z ML = 

} H . 2046 

I RAIL IMG FUGF 


XM = 

A . 6 3 6 0 

ZN = 

13. 1010 

TRAILING SPIRAL COORD I NAT FS 

XM i\j = 

( . 9*2 3 

Z M N = 

1 6 . 6 f 2 3 

XI ( 1. )= 104 . 09^4 

X I ( M ) = 

f l .0? V* 




*** VHR T I CAL R t ACT [ I in s :- 





WK PL = 

-20 . 6 366 

wKr»L = 

4. 1 4 IK 

WKSL= -30. 1 61 6 

T n I a i. 

L FAD INI, - 

-64.623 3 

WKPT = 

. 6006 

W KG I = 

l 1 . 16H3 

WKS T= —SI ,26 H4 

TUTaL 

1 KA I L I N I j = 

-30.6U04 



SASHOK I 71 IMTAL 

R h AC r 1 1 INSSS 





HK Pl_ = 

26.6201 

HK L = 

-M . 1030 

TIITaL lfahimi. = 

1 / . 4 2 L l 

HK.PT = 

-13,1707 

HKT = 

-30.2603 

tiiTal Trail iml = 

-62.4211 


MfO. = -74.3244 MSL = 7HH.H741 MP I. =-302.6473 TOTAL L FA I) I MO = -34.1444 
MOT = 74.0H74 MSI = 390,7706 MPT = 194.7001 TIITaL ( KAIL 140 ' 420.4084 


TOTAL TOROlIF = 5H1.4427 LA-INCHES = HO. 3902 < O-MFT E RS 

C-M.FP = . 44 1 7 .7773 
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Table 1. 3. Horizontal test for - 102 and - 0. 75 (upper bound) 

3**1 MP l IT r> ftitUMp r Ft* S*** 

■j , p , x I t f* ) . sk = 94. noon i s.nnnn 101.9994 ./son 

alpha ,M'Wn, ,phi f c = n.noon .hsh* 41.1999 ,o/>ou 

« , n = 17. O'jf'n 17 . oouu 


***!> f RFSS J AIAMFT 


T mf T A ( LHP ) 

= 1 ?7 

.HHRP 

THF r A ( LL P ) 

= HP . 301 9 

1 H F r A ( M p ) = 

A P , 3 A ^ 0 

THFTMM > 

- H4 

.1430 

T H \- T A ( MM ) 

= S P.9 346 

I HP T A 1 MM ) = 

17 .SS47 

pn - 

PUPS 

P L i. = 

.8960 

HML ~ 1 

.07 HO 



PM = \ . 

07 9 1 

PMM = 

.887 \ 

PMM ~ 

-HA by 



NORMAL S T R P S S F $ 

AT L » M , 

, M = 

.P 3 7 < 1 . 

P 3P 3 

-A AMU 


T AM(iHNl 1 A). 

STRH.3SFS AT 

L ♦ M T M = 

. i s i 3 

7 01 V 

. 4P 16 



V<**G*tIJM*tTK I C u Ai< AMF TF4 S*** 


CFMIHR n- IMS r AMTAMPOUS RIKA ITMM 

XR = 

o. unuu 

7 H = 

J ( > . 4 6 P S 

TRAILING SPIRAL PHI> 

X P I = 

-p. v h h / 

/PK 

7 . 7 3 H 8 

L FADING SPIRAL poll 

XPP- 

-4. K>33 

7PP = 

11 .94 *4 

HI FI IRC AT I UN PIIIMT 

XM = 

-P. 90uj 

/ M = 

1 3. 44 HI 

L FADING FlKih 

XI. = 

~4. S9HU 

7 L = 

) 3 . ] 704 

L fading spiral gudko in at fs- — 

X^L = 

-4. 1 1 3^ 

7Mi. = 

16 . 7 SS 6 

trailing fih;f 

XM = 

] . H3UH 

/ N = 

}3 . 8 S 4 3 

TRAILING SPIRAL CHUR U I MA T FS 

XM !M = 

3. 4 S6 1 

/Mr.! = 

1 6 . P 7 1 P 

X I ( L ) = lOM.PASf) XI ( !M > = H3 

. P *74 




VERTICAL R F A C T [ 






w< P i.“ 

-P3 .8804 

9KGI_ = 

P.4838 

WKSI_= -27.2 74 1 

Till A i. 

L L A 0 I im 0 * 

-4 H . h / P 7 

WKPT = 

-3 .S70 1 

M KG T = 

S .OS fH 

WKS T= —40. Hi 1 1 9H 

HIT AL 

IKMI. )M0 = 

** 4 S .4 LP4 



**<‘Hi)KI7riNTfl.'L 

RFAO T ( 1 IN*** 





HK PL = 

14 * M/S 1 b 

HK L = 

- U -0S3S 

Till A 1 . LtAI) I Ml, = 

5 .9U30 

HK PT - 

-18.34 OH 

HKT = 

-PP. SH^4 

1 f IT AL 1 8 a 1 L 1 wo- 

-40.9030 


***MIMPimTS*»* 

MOL = -15.7603 MSL = 271. 1944 mPi. = -16i),H 970 TillrtL Fao[mO = 105.0470 
MOT = 7. 17HO MST = 27 7.47 80 w FT = 748.1724 TOTAL TkaIlInO* H43.H289 


TOTAL THROIIh = 448.8710 LS-pMCHFS = 84.704H KO-MFl'FKS 

FM.ff = .6897 .172 7 
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Table 1. 3a. Horizontal test for = 102 and = 0. 75 (lower bound) 

* [ :\ f Pill' A H A w l- l l- >i s o' 'r * 

i.i,P,xi = o^.noon « hum » luL.sss^-/ , hiiii 

AI.ihhv »<"» Am-i i \ t >j.u , r; = n.nnon 4-s.lwsv . 04oo 

P 7 t M ~ . I's.o^no 1 ?,ihjou 


* s;: rn y s S s ak a m r f i : K S - o' r 

I H)*Ta { L I : ^ ) = Mii.mIi 1 )] fHM‘«(LLH)= 1 hMAI«P)= 4/-. 3430 

|MI.-7/W'*U -IP, I 4 -ID TH^rfAf JMN) = I Ml- ( a ( roll ) ~ 3 3 . <'+'74 7 

<•' i i = ,^i)Mi P l L = ,/i/t'i 1 P hi I.. = . o ^ / h 

Pl« .0 ?*,_■> Pl*'M - . Mil'll = , ?4 3U 

•-ill p <i/i| Sf-A'SSfS aT I. r >-i P ■'■ = .14- Ms 1 *0*3/4 . I. OMs 

) A N«(i“ : 'i I T A i. ^ T K - S S • : s jv r I . , l '-'i » <\| ■= .')/?/ . h04A * ] A 4 M 


R It'll- TK } C PAKA-'-ih f hKS ! ^'v 


i y m T p ^ ( 

‘<4 1 

N S i A i\l T A >-i P ( ) 

ilS Kit 1 ATI (IN 

A A - 

0. 1 M)Lm> 

/h = 

1 0 . 

Ah/M 

TMM'.l »'i r. 

; s p 

[HAL H 1 1 L M - 

-- 

XM1 = 

-? . 7 

/HI- 

7 * 

/ MM IS 

1. 1 fl!l|5lf, 

SP [ 

KAL MI'U.1 -- 

- 

xpV = 

-4. 

/H^ = 

1 1 - 

SA / 4 

\< i i : h i ; a ( 

i 1 ])■! 

H>: ) 1 Ml 


X * v - = 

. SOL ^ 

/M = 

1 M • 

IMS 1 

LF-Ao I 

1 1 c; 

t 


XL = 

— 4 • H A 4 

/I - 

1L 

0 I'? A 

map f m i ; 

Slj i 

HAL 0 1 »> i i-i ' J l ■ 

■J A T MS--*- 

X 1 **L - 

-a. 3,Vu 

/ M L = : 

1 l * 

v/HM 

( HA 1 i. ] A,r: 


i,i- 


X “i = 

S . S / Uh 

/M - 

1 M , 

H/?V 

1 HA 1 i. 1 N ( ■ 

. 3P 

i Hal (.■' ’i > vn 

| MA f H s 

X'M = 

/, 4 S 

/ i.i iM - 

1 A . 

] 04? 

* T ( L > “ 

1 10 

. 4 / A M 

X I 1 M 5 = M. 

A A/ 1 






• v * VmMIC a I. ■< I- A ( ; 1 r i } ;\i -? =;- 


w < P |_ = -2.1 . 3s<, 3 .M K(; l_ = 4.715/ WKSL = -36. ‘31 15 Till A', 1.5/inlwi, = -52.513/2 

Vi<PT= -.4757 -iKIi f = M.466K wHSTs -52. 7743 rn (ill. fKAlefi' 1 !, ■ -41.0tM5 

is I 7 1 1 i' 1 1 A •. Sr/WIT ini'.)** 1 * 


H<r>l.= 73.3 344 1 1 K I. - -H.31L2 nil «L Lbfln [m, = 14.H232 

lK»Ts -I 4. lots 3 H5T = -3S. <>2/0 llllAl. r Kf, 1 1. I rlt, - -44.H732 


- ti-ihM r S=:= 

1 *1. - -2s. >rnl 3 -'SI. = 250.5162 MW. =-265.0474 nil At '.rttOl'Ml, = 5.4404 

b(;l = 31.7610 ..ISI = 364 . 733 / H-> I . 210.1747 Till A I. l5AIU'H,= 445.6100 


(MAI. f. 12011.. : 6(14.1004 Ln- 1 I'itlH 64 = 43.6701 KO-'ll- l> KS 

1,6 0 - ,01.62 .3/2 i 
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Table 1.4. Horizontal test for £; ^ = 104 and “ 0. 75 (upper bound) 

***IMPIIT PARAMETERS*** 

W.P.XI <MI,SK = 04.0000 35.0000 103.0000 .7500 

ALPHA .GAMMA , PH! ,C = 0.0000 .05H4 34.0000 .0400 

R.R = 13.0500 IK. 0000 


***STRESS PARAMETERS*** 

THETA ( LOP ) = 147.0000 lHEfA(Ll.P)= 92.4864 1HETAIMP)= 47.6787 

THF T A (Ml = 102.6787 THETAINN1 = 38.3406 f HE T A (NO I = 14.7468 

PO = .2009 PLL = .7616 PMI. = .9852 

PM = .9567 PMN = .6260 RNM = .3577 

NORMAL ST RES St 5 AT L.M.N = .2845 1.1647 .1637 

TANGENTIAL STRESSES AT L.M,N= -.1940 .4635 .1742 


***OEOMETRIC PARAMETERS*** 


center of instantaneous rotation — 

XR - 

o. ooou 

78 = 

10.4625 

TRAILING SPIRAL POLE 

XP 1 = 

-2. 1518 

ZPl = 

H .0994 

LEADING SPIRAL POLE 

XP2 = 

-5.4226 

2P2 = 

11 . 2H67 

BIFURCATION POINT 

XM = 

-3.3748 

ZM = 

13.5356 

LEADING EDGE 

XL = 

-5.4894 

ZL = 

J 2 . H24*3 

leading spiral coordinates — 

xml= 

-5. 65Ub 

ZMI_ = 

16.5409 

TRAILING EDGE 

XM = 

4.359b 

ZN = 

13.2512 

TRAILING SPIRAL COORDINATES 

XMN = 

7.4397 

ZMN = 

1 6.6879 


XI ( L 1 = 113.1730 XI ( N ) = 71. 78 80 

«>#;s VERT I CAL REACTION*** 


WKPL= -17.3331 WKGL = 3.4031 WKSL= -31.6086 TOlfll LEADING = -45.5386 

WKPT= 2.2742 WKGT = 13.0597 WKST = -63.6576 TOTAL TRAILING* -48.3235 

***HOR I 70 NT A L REACTION*** 


HK PL* 23.1116 HKL = -7 .7770 TOTAL LEADING = 15.3346 

HKPT = -16.68HL HKT = -34.6464 TOTAL TRAILING* -50.3346 


***MUMENTS«*» 

MGL = -22.3480 MS L = 272.3250 MPl =-244.5668 TOTAL LEADING = 5.4101 

MGT = 32.6770 MST = 332,9537 MPT = 242,4732 TOTAL TRAIL I NG= 608.1040 


TOTAL TORQUE = 613.5141 LR-INC.HES = B4.82L5 KG-MFTERS 

FM.ET , THRUST EFFICIENCY: .6238 .3723 .5968 
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Table 1.4a. Horizontal te st for = 104 and = 0.75 (lower bound) 

#*#INPWT PARAMETERS*** 

W,P,XIIM),SK = 94.0000 35.0000 103.9999 .7500 

ALPHA, GAMMA, PHI ,C = 0.0000 .0584 34.9999 .0000 

R.n = 13.9500 12.0000 


***STRESS PARAMETERS*** 

THETA ( l_np )= 14H.4617 THE T A ( LL P ) = 95. 3393 1 HE T A [ HP 1 - 47.6787 

THETA(M) = 102.6787 THETA ( NN I = 32.7638 THETA (NO 1 - 21.2900 

PO = .2009 PLL = .7183 PMl = .9691 

PM = .9313 PMN = .5446 PMM = .2659 

NORMAL STRESSES AT L.M.N = .2444 1.1315 .0933 

TANGENTIAL STRESSES AT L,M,N= -. 1379 .4512 . 1253 


***GEOMF TRIG PARAMETERS*** 


CFMTER OF INSTANTANEOUS ROTATION 

XR = 

0 . oouo 

ZR = 

10.4625 

TRAILING SPIRAL POLE 

XP1 = 

-2.1518 

ZPl = 

R .0994 

LEADING SPIRAL POLE 

X P2 = 

-5.4226 

IP 2- 

1 1 . 2867 

BIFURCATION POINT 

XM = 

-3.3748 

ZM = 

13.6356 

LEADING EDGE 

XL = 

-5. 5887 

2.1 = 

12.7815 

LEADING SPIRAL COORDINATES 

XM l = 

-6.0313 

Z ML = 

16.7657 

TRAILING EDGE 

XM = 

6.3083 

ZM = 

12.9005 

TRAILING SPIRAL COORDINATES — 

XMN = 

8.8594 

ZMN = 

15.1858 


X1IL)- 113.6174 XT f M > = 6 7.6336 

*** VERTICAL REACTION*** 


O K P L = -16.4021 WKGL = 4.0320 WKSL= -34.5012 TOTAL LEADING = -46.8713 

NKPT= 2.7157 NKGT = 16.0181 WKST = -65.5958 TOTAL TRAILING* -46.8619 

* * * H 0 R I 7 O N T A L REACTION*** 


HK P L = 25.0784 HKL = -6.9340 TOTAL LEADING = 18.1443 

HK PT = -13.2243 HKT = -39.9200 TOTAL TRAILING* -53.1443 


***MOMENTS*** 

MOL = -25.5632 MSL * 280.4473 MPL =-277.4984 TOTAL LEADING = -22.6143 
MGT = 62.2253 MST = 364.2121 MPT = 207.2344 TOTAL TRAILING* 623.6719 


TOTAL TORQUE = 601.0575 LB-INCHES = R3.0993 KG-MFTERS 

EM, FT, THRUST E FE I C I ENC Y = .6111 .3723 .6092 
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Table 2.1. Wheel performance on slope for = 110 and = 0.55 (upper bound) 

***INPUT PARAMETERS*** 

w.p.xiimi.sk = ion. 30 oo o.nooo 109.9999 .5500 

ALPHA , GAMMA, PHI ,C = 20.4159 .0564 43.1999 .0600 

R ,B = 1 3. 0500 12.0000 


***STRESS PARaMFT ERS*** 

THF TA ( LOP I = 135.4904 THE TA { LL P I = 91.8329 THETA ( M P ) = 61.6863 

THETA(M) = 10H.48R3 THETAINN) = 81.2609 THETA INf) I = 43.3603 

PO = . 2025 PLI. = . 64 75 PMI. = 1.0555 

PM = 1.0314 PMN * .9135 PNN * .1011 

NORMAL STRESSES AT l_«M,N = .4372 .95 94 .5715 

TANGENTIAL STRESSES AT L»M f N= -.4654 .706U .4754 


***GEOME TR | C PARAMETERS*** 


CFNTFR OF INSTANTANEOUS ROTATION 

XB = 

-2.6764 

ZB = 

7. 1905 

TRAILING SPIRAL POLE 

XP1 = 

-5.6477 

Z P 1 = 

1.2107 

LEADING SPIRAL POLE 

XP2 = 

-10.374 i 

ZP2* 

6. 14 I 7 

BIFURCATION POINT 

XM = 

- 4 * 0442 

ZM * 

10 .6209 

LEADING EDGE 

XL = 

-10.4L61 

ZL = 

9.2793 

LEAniNG SPIRAL COORD] NATES 

XML = 

- 10. 5? 74 

ZMI. * 

1 2. 7547 

trailing edge- 

XM = 

-4.0304 

ZN = 

13.3549 

trailing spiral cooroinates 

XMN= 

-3.5423 

ZMN* 

16.5337 


XI (LI* 117.6875 XI<N>= 86.3796 

*** VHR T ! CAL REACTION*** 


WKPL* -19.2240 WKGL - 1.9010 WKSL* -21.9923 TOTAL LEAPING = -39.3153 

WK PT = -13.2269 WKGT - 6.1382 WKST* -54.3763 TOTAL T RA I L 1 NG = -61 .4650 

***HORIZONTAL REACTION*** 


HRPL* 19.0691 HKL = -3.2536 TOTAL LEADING * 15.6154 

HKPT* -16.2901 HKT = .4747 TOTAL TRAILING* -16.6154 


***MOMENT$**« 

MGL = -21.1603 MSI * 261.8661 MPC = -9.9469 TOTAL LEAPING = 230.7766 
MGT = -38.5212 MST = 244.5758 MPT * 294.5266 TOTAL TRAILING* 500.5603 


TOTAL TOROIIE = 731.357? LH-INCHES = 101.1139 KG-MET6RS 

EM, FT. THRUST EFFICIENCY* 1.0140 .3722 .3670 
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Table Z. la. Wheel performance on slope for = 110 and Sj_ = 0.55 {lower bound) 

«**inpdt p apame tek 


w,p,xi(M),sk = mo . 3000 o.oono 109.9999 . 5300 

ALPHA .GAMMA .PHI , C = 20.4159 .0684 43.1999 .0000 

R.H =* 13.9*100 12.0000 


***Sf8ESS oAPAMFTEKS*** 

THETA! LC1P 1 = 137.2*309 THETA|LLP)= 100.7263 IHETA(MP)= 61.6833 

THFTA(M) = 108. 48 H 3 fHETA(NN) = 78.0667 f HE TA (NTT ) = 6 1.7660 


P() = 

. 20 25 

P IL 

- .6706 

PML = 

.9648 


PM = 

.9181 

PMN 

= .7337 

PNN = 

.4207 


NORMAL 

STKfSSPS 

AT I 

,M,N 

.24 19 

. 84 70 

• 317M 

tangential stresses 

AT 1. 1 Mi N = 

-.2787 

.6284 

.2853 


***f;FLlMFTRlC PAR AMP TbRS*** 
CFNTFR OF INSTANT ANFDUS ROTATION XR = 

- 2 . 

78 = 

7. 1905 

TRAILING SPIRAL PHL F 

XP 1 = 

-b.8977 

7P1 = 

1 .2107 

LFAOINf; SPIRAL ROLF 

X P^ = 

-10.3797 

2P2* 

8. 14 17 

H I FI IRC AT I ON POINT 

XM = 

-9.044*! 

2M = 

10.6209 

L FA 0 I NG FDGt 

XL = 

-10.S634 

7l = 

9.1L13 

L F AO 1 NG SPIRAL COORDINATFS 

XM L = 

-11.373b 

7 ML = 

13.3880 

TRAILING FOGF 

XN = 

-2. 31 

7N = 

13.7564 

TRAILING SPIRAL COORUINATFS 

XMN = 

- 1 . 10b 9 

ZMN = 

17.9838 


XI I L > = 1)8.8061 XIINI= 79.1403 

*4** VERTICAL H FACT t ON* 44- 


WK PL = -17.3385 WKGL = 3.1863 WKSL = -27.6860 TOTAL LEADING = -41.7393 

W K. P T = -11.0964 WKGT = 1 1.0070 WKST= -68.5678 TOTAL T 8 A I L 1 NG = -58,6663 

***HI)Kr7tlNTAL REACTf l)N»*** 


HK PL- 22 .9886 HKL = -1.2492 TOTAL LEADING = 2 1.7393 

HKPT= -16.4917 HKT = -6.2476 TOTAL TRAILING* -21.7393 


Y^MOMENT S*** 

MC,L * -34.0184 MSL = 307.4669 MPL = -70.9631 TOTAL LEADING = 202.4863 
MGT = -66.1 329 MS T * 242. 1864 MPT = 283.2736 TOJAL THA 1L ING* 469.3271 


TOTAL TORDlie = 6 7 1.8124 L8-INCHES = 92.8815 KG-METEKS 

FM.ET, THRUST EFFICIENCY .9315 .3722 .3995 
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Table 2.2. Wheel performance on slope for = 110 and - 0.65 (upper bound) 

***INPIJT PARAMETERS*** 

W,P,X|(M),SK = 100.1000 0.0000 109.9999 .6500 

ALPHA, GAMMA, PHI ,C = 20 .ALSO .0586 63.1999 .0000 

R,8 = 13. 9500 1H. 0000 


***S TRESS PARAMETERS*** 

T HFT A { LDP > = 136.7667 THE f A ( UL R ) = 86.0939 THE T A ( MP 1 = 70.1512 

THFTAIMl = 110.0512 T MET A { NN I = 88.2381 1HE1A(N0)= 46.8502 

PO = . 2025 P LL = 1.0312 PM 1 . = 1.1467 

PM = 1.1313 PMN = 1.0572 RNN = .8398 

NORMAL STRESSES AT L,M,N = .<5501 1.285U .7561 

TANGENTIAL STRESSES AT -.6263 .7637 .676H 


* A * Gfc (TM E T R 1 C P AR A M F T FH S* 6 * 


CENTER OF INSTANTANEOUS ROTATION 

XB = 

-3* 1630 

ZB - 

H .6979 

TRAILING SPIRAL PULF 

XP 1 = 

-5.1566 

ZP1 = 

2.9751 

LEADING SPIRAL PULE 

XP2 = 

-ll).?565> 

ZP? = 

7.2615 

BIFURCATIUN POINT 

XM = 

-9.044^ 

ZM = 

10.6209 

LEANING EDGE 

XL = 

-10.0497 

ZL = 

9.6769 

L FADING SPIRAL COORDINATES 

XM L = 

-9.8667 

zml = 

11.8076 

trailing edge — 

XN = 

-4.8458 

ZM = 

13.0812 

TRAILING SPIRAL CnORDlNATFS 

XMN = 

-4,7345 

ZMN = 

1 6.7006 


X T { L ) = 115.6725 XI(N)= 39.0106 

***VERTIGAL REACT ( (INS'"?* 


WKPI_ = -15.0300 WKG!_ = .7116 WKSL= -13.6206 TOTAL L FADING = -27.7392 

WK PT- -19.9701 WKGT = 6.3686 WKST= -59.17(10 TTIlAL T KA I L I NG = -72.7717 

***H0R17.0NTAL REACT! ON* ** 


HKPL= 11.9736 HKL = -Z.0915 TOTAL LEADING = 9.8821 

HKPT = -19.7695 HKT = 9.8676 TOTAL TRAILING^ -9.8821 


***MDMENT 5**“ 

MGL = -8.5330 MSI = 157.1366 MPL = 20.8128 TOTAL LFADING = 169.6166 
MGT = -62.8072 MST = 216.5360 MPT = 391.0730 TOTAL TRAIL ING= 662.8018 


TOTAL TORQUE = 732.2163 LH-INCHF.S * 101.2327 KG-METERS 

FM.Ff .THRUST E FF I C I FNCY = .8690 .3722 .6332 
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Table 2. 2a. Wheel performance on slope for = 110 and = 0.65 (lower bound) 

*®*”IMP(jT PARAMETERS*** 

W.P.XIIMl.SK = 100.3000 0.0000 100.4900 .6500 

ALPHA.GAMMA.PHI ,C * 20.4164 .0584 A3. 1444 .0600 

R.H = 13.4500 12.0000 


***STRESS PARAMETERS*** 

1 El E T A ( L 0 P )= 134.5818 THET A( LLP 1 = 89.0328 TFIETA(MP> = 70.1512 

THE TA < M I = 116.9512 EHEfAlNNI = FA. 9370 THE TA ( NO 1 - 48.1827 

PI) = .2025 PEL = . 9 0 1 A PML * 1 .0527 

PM = 1.0320 PMN = .9366 PIMM = .6766 

NORMAL STRESSES AT L,M,N = .A473 1.1661 .6929 

TANOFNTIA'. STRESSES AT L.M.N* -.5148 .6 782 -A6?l 


***GEUMETR I C PARAMETERS*** 


CFNTHR OF I N ST AN TANFUl IS ROTATION 

XB = 

-3.1 630 

28 * 

8.4979 

TRAILING SPIRAL POLL 

X P 1 = 

-5.1566 

ZP1* 

2.975 1 

l fading spiral POLL — 

XP2 = 

- 1 0.2669 

2P7 = 

l. 2615 

bifurcation point — 

XM = 

-9.0442 

ZM = 

10.6209 

L FADING PDG! 

XL = 

-1 0.2191 

ZL = 

9.4958 

lfaoing Spiral coordinates — 

XML = 

-10.1747 

ZML * 

1 2 .1278 

TRAILING EDGE 

XN = 

-4.2429 

ZN = 

13.2891 

TRAILING SPIRAL COORD [NATES 

XMN = 

-3.8774 

ZMN* 

17.4139 


XI ( L ) — 116.6863 XI(N)* 87.2911 

*#*V t R T I CAE REACTION*** 


WK PL = -15.6411 WKGL* 1.1163 WKSL* -1 6.0262 TOTAL LEADING = -30.5610 

WKPT* -18.31A1 WKC, T = 8.275A WKST = -59.6836 TOTAL TRAILING* -69.722 3 

***HDR I ZONTAL REACTION*** 


HK PL = IA.1199 HKL * -1.7898 TOTAL LEADING = 12.3300 

HKPT = -19.8063 HKT = 7.A762 TOTAL TRAILING* -12.3300 


***MOMFNT S*** 

MGL * -12.8101 MSL = 183. 1972 M PL = 6.3AA5 TOTAL LEADING = 176.7316 

MGT = -52.2301 MST = 212.82AA MPT = 380,4392 TOTAL TRAILING* 541.0335 


TOTAL TORQUE = 717.7651 LB-INChES = 

FM, FT, THRUST EFFICIENCY* .8421 


99.2347 KG-METERS 


,3722 .4420 
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Table 2.3. Wheel performance on slope for = 105 and = 0.65 (upper bound) 

*«*INP(/T PARAMETERS*** 

W,P,XI <M) , SK = ion. 3000 0.0000 104.9999 .6600 

ALPHA, GAMMA, PH! ,c = 20.4164 .0684 43.1449 .0600 

K , R - 13. 9600 12.0000 


***S1RESS PAR AM F T EKS *** 

TH1 TA ( i.flP ) = 140.1,762 THtzlAlLL P)= 96.7496 THETAIMP)* 69.7416 

THETAIM) = 106.6416 T H t r A ( NN ] = 76.6602 IHflMNIIP 37.2337 

Pll = . 2025 PLL ' . 8 64 7 PML = 1.0901 

PM = 1.066? PMN = .9376 PNN = .7373 

NORMAL STRESSES AT L,M,n = .299? 1.0610 .6678 

TANGENTIAL STRESSES AT L,M,N- -.3172 ,714b .4913 


***GEUMETR1C PARAMETERS*** 


CENTER UF InSTANTANEUUS RfJT A T I UN 

XB = 

-3. 1630 

L B = 

8 .4979 

TRAILING SPIRAL POLE. 

XP 1 = 

-5. 8689 

Z Vi = 

3.8767 

LEADING SPIRAL POLL 

X P^ = 

-9. 3819 

ZPZ* 

9. 1441 

BIFURCATION POINT 

XM = 

-8.0841 

ZM - 

1 1 . 3687 

LEANING FOGC 

XL = 

-9.6078 

ZL * 

10.2080 

LEADING SPIRAL COORDINATES 

XM L = 

-9.9371 

ZML = 

13.8361 

TRAILING EDGE 

XN = 

-3.6743 

/N = 

1 3 .4829 

TRAILING SPIRAL COORDINATES 

XM N = 

-2.9136 

Z MN = 

1 6 . 288 3 


XI ( L 1 = 1 1 2.6601 XI I N ) = 84.4613 

*** VERT! CAL REACTION*** 


VtKPL* -18.7890 VOTGL* 2.3176 VOSS L* -26.9640 TOTAL LEADING = -43.4364 

W K P T = -1 1 .0238 HKGT* 4.9676 WKST* -61.2681 Tf)7 AL TRAILING* -67.3246 

***HHR 1 ZONTAL REACTION*** 


HRPL* 22.2113 HRc = -3.0039 TOTAL LEADING = 14.2073 

HKPT* -16.1693 HKT = -3.0380 TOTAL TRAILING* -19.2073 


***MDMENTS**» 

MGL = -22.9016 MSL = 288 .7918 MPL = -86.8012 TOTAL LEADING = 180.0890 
MG T = -27.3363 MSI = 268.0377 MR 1 = 27 7.2376 TOTAL TkAlLlNG* 607,9399 


TOTAL TllROUt = 68R.02R9 *.8- ! NCHR S * 96.1236 KG-MfcT E RS 

EM, FT, THRUST EFFICIENCY* .H072 .3722 .4611 
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Table 2.3a. Wheel performance on slope lor £ ^ = 105 and ~ 0.65 (lower bound) 
***INPUT PARAMETERS*** 

W,P,xr(M) f SK = 100.3000 0.0000 104.9999 .6500 

ALPHA .GAMMA , PHI ,C = 20.4159 .0584 43.1999 .0600 

R,B = 13.9500 12.0000 


***STRES$ PARAMETERS*#* 


THETA ( LOP 1 = 144.6045 

THETA(M) * 106.5416 

THETA! LL P ) = 111.9612 ThFTAIMP!* 

THETA(NN) = 60.8132 THETA! NO ) * 

59.7416 

56.5732 

PU = .2025 PLL = 

PM = .8925 PMN = 

. 5905 PML = .97 42 

.6214 PNN = .2 327 



NORMAL STRESSES AT L » M» 
TANGENTIAL STRESSES AT 

N = .1238 .8866 

L t M » N= - • 0 344 • 6 082 

.1220 

.1522 



***GEOMETRIC PARAMETERS*** 
CENTER OF INSTANTANEOUS RUTATIUN X8 = 

-3. 1630 

Zb = 

8.4979 

TRAILING SPIRAL ROLE 

X PI = 

-5.8589 

ZP1 = 

3.8767 

LEADING SPIRAL ROLE 

X P2 = 

-9.3819 

ZP2 = 

9.1^41 

BIFURCATION POINT 

XM = 

-8.0841 

ZM = 

1 1 . 3687 

LEADING EDGE 

XL = 

-9. 7268 

IX = 

9 . 9995 

LEADING SPIRAL COORDINATES 

XM L = - 

1 1 .6486 

Z ML = 

14.7655 

TRAILING EDGE 

XN - 

-.2333 

ZN = 

13.9480 

TRAILING SPIRAL COORDINATES 

XM N = 

2.2053 

ZMN = 

18.3139 


XI(L)= 113.7920 XI(MI= 7(1.5424 

*** VERTICAL REACTION*#* 


WKPL= -14.2141 WKGL* 4.9177 WKSL* -36.9407 T(!TAL (.FADING = -46.2371 

WKPT* -6.3939 WKGT* 14.5754 WKST* -62.3644 TOTAL TRAILING* -54.1829 

***H0 R I ZON TAL REACTION*## 


HK PL = 29.1273 HK.L = 2.5820 TOTAL LEADING = 31.7094 

HKPT* -14.7801 HKT = -16.9292 TOTAL TRAILING* -31.7094 


***MDMENTS*** 

MGL = -44.7884 MSL = 338.9945 MPL =-213.3103 TOTAL LEADING = 80.8957 

MGT = -45.3722 MST = 317.9378 MPT = 235.8248 TOTAL TRAILING* 508.3904 


TOTAL TORQUE = 589. 2861 Lrt- 1 NCH ES = R1.47LB KG-METEKS 

FM.ET, THRUST EFFICIENCY* .6913 .3722 .5383 
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Table 2.4. Wheel performance on slope for = 105 and s k = 0. 70 (upper bound) 

##*input parameters*** 

W ,P ,X1 I M) ,SK- = 100.3000 0.0000 104.9999 .7000 

ALP HA .GAMMA, PHI, C = 20.4159 .0584 43.199V .0500 

R » H = 13,9500 12.0000 


***3TRESS PARAMETERS*** 

THFTAtLOPt* 140.1015 THETA! LLPI* 95.0272 THETA (MR ) = 54.6400 

THE T A ( M I = 111.4400 THE T A ( NN 1 = 79.0130 THETA I Nfl ) * 39.6779 

PO * .2025 PLL * .8875 PML = 1.0926 

PM * 1.0613 P MIX = .956? PNN = .7353 

NORMAL STRESSES AT l.M.N = .3302 1.1885 .5801 

TANGENTIAL STRESSES AT L.M.N* -.3545 .7009 .4950 


***GEUMETR I L PARAMETERS*** 


CENTER OF INSTANTANEOUS ROTATION 

XB = 

-3.4063 

2B = 

9.1516 

TRAILING SPIRAL POLE 

XP1 = 

-5.4884 

ZP1* 

4.75H8 

LEAPING SPIRAL POLE 

XP2 = 

-9.3246 

ZP2 = 

8.7518 

8 IEORC AT I ON POINT 

XM = 

-8.0841 

ZM = 

11 . 3687 

leading foge — 

XL = 

-9.4567 

ZL = 

10.2553 

LEAPING SPIRAL COORDINATES 

x Ml = 

-9. 7420 

ZML = 

13.4988 

TRAILING E OGF 

XN = 

-3.8068 

ZN = 

13.4205 

TRAILING SPIRAL COORDINATES 


-3.1b5b 

ZMN* 

16.6198 


X I ( L 1 = 112.2640 X I < M ) = 85.4203 

#**VERTICAL REACTION*#* 


WKPL= -17.5474 WKGL= 1.8740 WKSL* -23.7603 TOTAL LEADING = -39.5337 

WKPT = -13.3298 WKGT = 5.6038 WKST = -53.061b TOTAL TRAILING* -60.7876 

***H0R1 70NTAL REACTION*#* 


HK PL * 19.6827 HKL * -2.2556 TOl At LEANING * 17.4271 

HKPT* -17.9479 HKT = .5208 TOTAL TRAILING* -17.4271 


##* MOMENT S#*# 

MGL = -18.6533 MSL = 248.9197 MPL = -65.3835 TOTAL LEAPING = 164.8829 
MGT = -31.2311 MST = 237.8679 MPT = 317.2553 TOTAL TRAILING* 523.8920 


TOTAL TORQUE = 688.7750 L8- I NCHES * 95.2267 KG-METERS 

FM, FT, THRUST EEEICIENCV* .7503 .3722 .4960 
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Table L. la. Wheel performance on slope for | ^ - 105 and = 0. 70 (lower bound) 

**•» INPUT PARAMETERS*** 

W,P,XI(M).bK = 100.3000 0.0000 104.9994 .7000 

AL e EIA, GAMMA, PHI ,C = 20.4159 .0544 43.1999 .0400 

R.R = 13.9500 12.0000 


**#stress parameters*** 

THFTAILOPl* 142.3620 THHTA(LLP)* 106.4124 THETA IMP ) = 64.6400 

THFTA(M) = 111.4400 THErA(NN) = 61.2304 T HE TA 1 NO I = 50.6741 

PO = .2025 PLL = .6581 PMI. = .9738 

PM = .9 102 PMN * .7086 PNN = .3485 

NORMAL STRESSES AT L.M.N = .1604 l.OlOl .2184 

TANGENTIAL S IRESSES AT L,M,N= -.121 f .6011 .2292 


***OEOMETR 1C PARAMETERS*#* 


ONTPR nf I N S T AN TANFOUS « 0 T A 7 TUN 

XR = 

- 3 . 4 06 3 

ZH = 

9.1516 

TRAILING SPIRAL PCI Ur 

XPl - 

— 5 .48b*t 

ZPl * 

4.7588 

L h AD ! NG SPIRAL PULE 

XP2 = 

-9. 374*3 

ZP2 = 

8.7518 

H I F l IRC A T I ON PUINT 

XM 

-8.UB41 

ZM = 

11 .3687 

LEADING FHGfc 

XL = 

-9.699 / 

ZL = 

10.0258 

LFAHING SPIRAL CHOROJNATFS 

XML = 

-10*9471 

ZML = 

14.2607 

TRAILING FDGF 

XN = 

-1.6729 

ZN *= 

1 3.8493 

trailing spiral cuoroinates — 

XMlSj = 

.1837 

ZMN = 

1 8 . 2728 


X I ( L ) = 1 13 .6369 XI(N)= 76.471 8 

***VERT I CAl REACTION*** 


WKPL- -15. 1587 WKGL* 3.6518 WKSL* -31.0287 TOTAL L FAD I NO = -42.5 350 
WKPT = - 9.1669 WKC.T* 12.1788 WKST* -60.4015 TOTAL TRAILING* -57.3926 

* * *H0 R I 7 1 1 N T A L REACTION*** 


HKPL* 24.8834 HKL * 1.1439 TOTAL LEADING * 26.0273 

HKPT* -17.3956 HKT = -8.63L7 TOTAL TRAILING* -26.0273 


***Mf)MENT 5*** 

MGL = -33.9626 MSL = 292.4108 MRL =-148.4692 TOTAL LEADING * 109.9789 
MOT = -48.9456 MST = 2R5.0673 MPT = 289.6205 TOTAL TRAILING* 525.7422 


TOTAL TOROIlE = 635.7212 L8-INCHES = 87.891 7 KG-MFTERS 

FM,FT,THROST EFFICIENCY* .6925 .3722 .5374 
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Table 3.1. Wheel performance for Apollo LRV for £ - 105 and = 0.85 (upper bound) 

*** INPUT p ATAMFTFUS*** 

W,P .XT (M>, SK = 60.0000 0.0000 1 04 . OHHH .05 00 

ALPHA .GAMMA, PH] ,C = 0.0000 .0100 32 . 9HHH .0501) 

R.B = 16.0000 10. 01)00 


*** b T R t. S5 JAfAHHRS*** 

T HF T A ( l.flP ) -■ 15 2.0047 TH5rML'.PS= 107.737? 1HFT6(MP(= 46 . H V 2 2 

THFTMM) = 122.R722 ThFTAINN) = 64.4247 )HPTA(f|lll= 24.3704 

PO = . 17i 00 Pl.t = .4664 PM'. = .4254 

PM = .6057 PMN = .4771 HNM = .414? 

NI1KMAL S ' P-4',- s AT L.M.N = .1414 .6H63 .3046 

TANGENTIAL STRESSES AT L,M,N= -.0674 ,00V« .2241 


*»#G EUME T« I C PARAMETER S#«» 


Cf-MT^R Hi- IimST ANTANFIJUS RO f T !(tM- 

— XH = 

0 . OlHJU 

/R = 

} 3 .6000 

TRATLTNG SPIRAL PULP 

XP L = 

- I . 2 04 3 

/Pl = 

10 .OU)7 

L FA H f MG SPIRAL P ( l L F 

X P2- 

- 6 . A 2 | | 

/P2 = 

10 . 3 64 2 

3 I FI IRC AT I ON POINT 

XM = 

-A. 1411 

/ M e 

1 3 . 4 54 H 

LFAH I MG Fnr,F 

XL = 

- 7 . *> 7 7 f 

7 L - 

\4 .001 f 

L F A D I N G SPIRAL Cl KIR DINATFS 

XM L = 

-9, 062 l 

l Mf_ = 

] 8.8 <9U 

TRAIL IMG Fl)(t*= 

XM = 

1 . 2092 

/ N = 

] 3 .9342 

TRAILING SPIRAL COOK!)! NAT PS 

XMN= 

3. 32 33 

/mm = 

20.3762 


XI ( L > = 1 1 B . ?6 67 XI ( N ) = H 6. 66 4 6 

##* VhRT I CAL REACH MN*«* 


VKPL = 
WKPT = 

— 6 .7 076 
-A. 0106 

WKGL = 1,1136 

WKGT= 1.7418 

WKSL = 
wKST = 

-21 . 9RU3 
-2 9.8 7h 3 

TOT A l. 
TOTAL 

LPAlMMG = 
TRAIL IM,= 

-2 7 . 66 2 3 
—32. 1433 



7UNTAL 

nmcn™«7? 




HK Pt. ~ 
HK PT = 

1 6 » 34 02 
~ ] 4 . 8678 

HK L = 1.8722 

HK T = -3* 3306 



total 
nn al 

lpaoing = 
1 RAIL 1 MG = 

18.2134 

-18.2184 



***MOMp 

NTS* * * 





MGL ~ 

MGT “ 

-7.2439 

2326 

MSI = 138.9629 
MST = 10.2989 

MPl = 
MPT = 

-21 3. 3682 
261. 60 l H 

nn al 
TOT al 

L FADING = 
TRAIL IMG = 

-61.8312 

271.6481 

TOTAL 

T GROUP = 

209.7969 LR- 

l NC h F S 

= 29. 

0033 KC 

,-MFThRS 


FM , FP 


= 

.237 1 

0. uooo 
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Table 3. la. Wheel performance for Apollo LRV for = 105 and = 0. 85 (lower bound) 

**#ENPHT PARAMFI FKS*** 

w,p, xhmi.sk = 60 . none) n.nono 104,4444 .hsou 

alpha, gamma, phi ,c = n.onoo .oirni 12.4444 .0600 

R,K = 16.0000 10,0000 


*** S T 3 F 3 S p A3 AM t- T HP S «* * 

T HF T A ( )_flP ) = 160. 8PM THl= r A ( LL P 1 = 122.4 616 1 FtF TA < MP ) = 66.8722 

THFTA(M) - 122.8722 TMfTMN«l = 64.2667 |IFTAI«(l)= 3.6.6041 

PI'l = .1600 P l_ L = ,4040 PML = ,484V 

PM = .460 3 PMN = . 3626 PNM 5 .3172 

NORMAL SIKFSSFS Ai L,M,N = .L103 .6014 .1474 

TAMOFMT I At. STRFSSFS AT l.,«,N= .0383 .0864 .1487 


* * # r, F (. | H f: T 8 I c O A i< A M F T F 3 S * >» 


CFMTPR HP I NS T AWT AKjPDU S RUUTIUN- 

— xr = 

o.ooou 

L R = 

1 \ ,/suuu 

trailing spiral phlf — • 

XP 1 = 

- 1 

/HI = 

10.4 1.07 

L FADING SPIRAL PULL 

XP^ = 

-h.4? 1 1 

/ V ?- 

1 0 • 

KIFUKCATIUN PUlNf 

X‘i = 

-a* 141 L 

i M = 

] h * A S4 H 

1 f-AMNG R()f»h 

XI. =: 

- h • '+ S b u 

L L - 

1-5 .*>7 12 

LEADING SPIRAL Cdl W ( > I M A H-S 

XA L - 

- 12. KJhV 

£tff. = 

L* . -W2? 

TRAILING Fl)Gh 

X N = 

3 • *3 S 6 to 

/M = 

1 S . SSKJS 

trailing spiral ci hjrij i ma t ls — 

XMlM = 

h.zshcj 

i Ri N = 

20 . 1 L b A- 


XI <1. ) = 121.4426 XI(NI= 7 7.0046 

* * * V K R T [ C A L RHACTION*** 


wkpi.= -3.0416 HK0L= 2.1307 WKSL= -27.4771 TUTaI. LFAOlMI. = -28.8380 

WKPF = .6744 WKGT = 3.6441 WKSI= -36.1432 ITITaL T MAILING^ -30.4144 

w*S--Ftf!K I 70NT4L R FACT 1 ON*** 


HK PL = 20.386 1 HKt = 7.5002 TOTAL LFAOIMO = 28.3864 

HKPT = -16.3401 HKT = -1 3.048? TOTAL T ka l !. I wij= -28.3864 


ttWWMUNFNTS*** 

M(1L = -10.3476 MSL = 164.3866 MPI. =-313. 7054 TlllAL L FADING = - 164.6474 
MOT = 7.8487 NST = 73.3065 HPT = 274.1733 TOTaL TkA!LInG= 360.0787 


TOTAL TOROtlh = 200 .3807 Lri-tNCHFS = 27. 7036 KG-MFTFRS 

F M , F P = .2468 0.0000 
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Table 3, 2. Wheel performance for Apollo I,RV for “ 105 and " 0. 90 (upper bound) 

***1NP0T PARAMFTPRS#** 

W,P,XI(M),SK ' AO .000(1 0.0000 104.4999 .4000 

ALPHA, GAMMA, PHI ,C = 0.0000 .0100 32.9994 .0500 

P. ,R = 15.0000 10.0000 


***STRFSS PARAMETERS*** 

THFT A ( LOP 1= 149.3241 THE T A ( U, P I = 97.2966 THFTA IMP )= 75.7046 

THFTA(M) = 1 32.7096 THE T A ( NN I = 79. 7621 T HE T A ( NO ) = 31.3747 

PO = .1690 PLL = .5499 PML = .5657 

PM = .5739 P MN = .5604 PNN = .5062 

NORMAL STRESSES AT L,M,N = .2061 .B094 .4919 

TANGENTIAL STRESSES AT L,M,N= -. 1360 ,00Mb .26 H4 


***GEUMF TR I C PARAMETERS*** 


CFNTFR OF INST ANT ANFUUS RH T AT I UN 

XB = 

o. ooou 

ZB = 

14 .4000 

TRAILING SPIRAL PIILF 

XP1 = 

-.6650 

ZPl = 

11.7107 

lfaoing spiral POLE 

x P2= 

-6.0214 

ZP2 = 

6.0725 

lUFHRCATjriN PUINT 

XM = 

-4. 141 1 

ZM = 

15.4546 

LEADING EDGE--- 

XL = 

-6.6394 

ZL = 

14 .4642 

LEADING SPIRAL CHORDI NATES 

XM L = 

-7.2572 

Z Mi_ = 

17.7234 

TRAILING EDGE 

XM = 

. 04 04 

ZN = 

15.9997 

TRAILING SPIRAL CHORD! NATES 

XMN= 

.4611 

ZMN = 

20 .6491 


XI (LI- 115.3066 XI ( M ) = 69.6762 

***VERTICAL R FACT I ON* ** 


WKPL= -7. 1043 WKGL= .4946 WKSL= -16.2660 TOTAL LEADING = -22.6960 

WKPT- -9.2000 WKF,T= 1.3934 WK8T? -29.5061 TOTAL TRAIL I NG = -37.3146 

***Hn K I 70NTAL RE AF.T I ON* ** 


HKPL= 11.5636 HKL = .7605 TOTAL LEADING = L2.3441 

HKPT= -16.5940 HKT = 4.2496 TnTAl TRAIL I NG= -12.344L 


***MF)MtNTS**« 

MGL = -3.6164 MSL = 120.1449 MRL =-136.5216 TOTAL LEADING = -19.9991 
MOT ? -1.7654 MST = -75.9649 MPT = 301.4636 TOTAL TRAILING?- 223. 7532 


TOTAL TORQUE = 203.7541 Lrt-INCHES = 76.1700 KG-METEKS 

FM,EP, = .2356 0.0000 


59 


JPL Technical Memorandum 33-477 



Table 3.2 a. Wheel performance for Apollo LRV for - 105 and - 0.90 (lower bound) 

***!NPIJT PARAMETERS*** 

WrPrXIlMlrSK = 60.0000 0.0000 Ml 4.9999 .VOW 

ALPHA , GAMMA, PHI , = 0.0000 .0100 32.9999 .(JSOU 

K,K = 10.0000 10.0000 


*»*i TRUSS PARAMETERS*** 


T HF TA ( LHP ) = 166 

THFTA < M ) = ] 32 

.1713 

.709b 

THt T A < LL F ) 
f H h T A ( NN ) 

= U9.7895 THSrA(MP) = 

= 46.7387 T HF T A 1 Nil ) = 

73.7096 

24.4482 

PU = ♦ 1690 

PM = .427 H 

PLL = 
PMN = 

. 38b 7 
.3715 

PMi. = .4642 

PNN = .27/0 



NORMAL SfRFSShS A T L * M i 
TANHFNT I A 1. STRtSSFS AT 

i N 

1. , M, N = 

.09 86 .5837 

.0010 .0064 

*170 1 
. 1 2 HO 



***GF|JMETRIC P ARAMFFERS*** 


center of ins r ANTANEOUS ROTATION — 

- XB = 

U . 00(H) 

/ 8 = 

14 .4000 

TRAILING SPIRAL PlILF 

XP 1 = 

-.6H3U 

Z H L - 

11.7 107 

L FADING spiral pole — 

XP2 = 

-6.0214 

ZP2 = 

(1.0723 

BI FIIRCAT TDN PIIINT 

XM = 

-4. 1411 

7 M = 

13 .4348 

LEADING EDGE 

XL = 

-a. 980 7 

l L = 

13.2416 

leading spiral ci ihr n i n at f s — 

XML = 

-12.2390 

/Mi, = 

I 8. 9683 

TRAIL! MG EDGE 

XN = 

3.12 71 

/N = 

13.6914 

TRAILING SPIRAL CflDR 1)1 NA T FS 

XM N = 

H . 7 06 9 

/MN * 

21.3168 


X I ( L ) = 

124 . 1454 

X I ( N 1 = 78. 77 90 






*** 

VERTICAL REACT I ONE 44 




WKPL- 
WKPT = 

-3.84 1 3 
-.1987 

WKGL- 
W KG T = 

2. 1766 WKSL= -26. I 374 
4.4923 WKS T = -36.3118 

TOTAL 

TOTAL 

L FADING = 
TRAIL I N G = 

-27.8020 
- 12. 0182 



***HPR I7DNTAL R E AG T I UN* ** 




HKPL = 
HK P T = 

19.0737 

-17.0687 

HKl. = 
HKT = 

8.3384 
-10. 3634 

total 

TOTaL 

L F A n I M ( > = 

TRAIL I N(i - 

27.4322 
-2 7 .4 422 




...MOMENTS*** 




MGL = 
MGT = 

-4.6617 

9.4906 

M S L = 
MST = 

166.3932 MPL =-267.8813 
24.2904 MPT = 31 8. 7603 

TOTAL 
1 0 T a L 

'.FADING = 
i RA I L IM(,- 

- 1 0 / .948 0 
3 6 2 . 6 4 1 2 


TOTAL THRONE = 244.393? I.H-tMCHFS 


33.H162 K G-ME 1 F R S 


FM , E P , 


. 7H30 0.0000 
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Table 4. 1. Wheel performance for Lunokhod-1 for - 112 and - 0. 80 (upper bound) 

***INPllT PARAMETERS*** 

W,P,XIIM),SK = 35.0000 0.0000 111.9999 . H 000 

ALPHA, GAMMA, PHI ,C = 0.0000 .0100 -42.9404 .0900 

R,H = 10.0000 6.50UU 


***STRFSS PARAMFT ER S'-** 


THF-T A 1 LOP ) * 155.8434 
THFTAfM) = 123.2469 

THETA(LLP) 
THE r A ( NN ) 

= 107.6?83 TH^TAIMP)- 

= 68.6^72 lHPTA(Nin = 

7 ] .2469 
24.4682 

PIT = .1690 PLL = 

PM = .53^5 PMrJ = 

.6043 

.5127 

PML = .5484 

PNM = .45 94 



NORMAL STRESSES AT l,m, 
TANGENTIAL STRESSES AT 

, N = 

L ,M,N= 

. 1 439 .7224 

.1447 .1207 

. 3H43 
*2502 



***G FUMETR I C PARAMO r EKS »** 


CENTER OF INSTANTANEOUS ROTATION 

X* = 

o.ooou 

ZK = 

8 .0000 

TRAILING SPIRAL POLE 

XP1~ 

-.8269 

Z P 1 = 

6.8672 

LEADING SPIRAL POLE 

XP2 = 

-6.2219 

ZP 2 = 

4.9247 

8 IHJRCAT ION POINT 

xm = 

- 3 . 74 6 0 

/M - 

9.2/18 

LEADING EnGfc 

XL = 

-6. 1601 

ZL = 

7.8773 

LEADING SPIRAL COORDINATES 

X Ml. s 

-7, 32 LY 

ZML = 

IL .6329 

TRAILING EDGE 

XM = 

.8982 

ZN = 

9.9696 

TRAILING SPIRAL COORDINATES 

XMN = 

2.6 64/ 

ZMN = 

14 . 20 34 


XI (U= 123.0258 XI 111)= 84.8463 

***VERT I CAL REACT ION*** 


WK PL = -3.4942 WKGL = .38 90 WKSL= -10.0841 TIDAL LEADING = -13.7781 

WKPT= —3.2561 W KG T = .9136 WKST= -18,4917 TOlAL TRAILING* -2 1 . 332 3 

***HUKl 20NT AL RLAC 1 MIN*** 


HKPL = P.7H46 HK L = 1.2993 TIITaL UP All 1 'ML = 10.01(39 

HKPT = -9.7117 HR 7 = -.3722 TOTAL TRAILING* -III. 0839 


***MOrtENTS*** 

MGL - -1.4460 MSL = 70.1045 M PI. = -61,8630 TIITaL LEADING = 6.2953 

MC,T = -.3179 MST = -31.3827 MPT = 112.0216 TOTAL TRAILING- 80.3208 


TOTAL TnROUE = 86 .6161 l.tt- INCHES = 1 1 .9761 KU-METEKS 

FM , FT , = . 3093 0 . OUOU 


JPL Technical Memorandum 33-477 



Table 4. la. Wheel performance for Lunokhod-1 for = 112 and = 0.80 (lower bound) 

«*<■ INPUT PARAMETERS*** 

W,P,XI(M),SK = 36 .0000 n .0000 111.9994 .8000 

ALPHA, GAMMA, PHI ,C " 0.0000 .0100 32.9949 .0500 

R » 8 = 10.0000 6.5000 


***STRFSS PARAMETERS*** 

THFTA I LOP ) = 188 .0952 THETAILLP)= 145. 0048 THFTA(MP)= 71.2469 

THFTA(M) = 128.2464 THfcTAINN) = 23. 3497 r HP i A ( NO ) = -2.5815 

PO = .1640 PLL = .4490 HML = .4347 

PM = .4812 PNN = .3895 POM = .3047 

NORMAL STRESSES AI L,M,N = . 1436 .6426 . 1376 

TANGENTIAL STR'SSFS AT L,M,N= .0873 .1087 .1343 


***GFUMFTRI C PARAMF T ERS *** 


CENTER nF INSTANTANFOOS ROTATION 

XK = 

o.oooo 

ZB = 

H .0000 

TRAILING SPIRAL POLF 

X P i a 

-.8254 

ZP1 = 

5,5672 

LEADING SPIRAL POLE — 

XP2 = 

-5.221V 

ZPA = 

4-. 924 7 

BIFURCATION POINT 

XM = 

-3.7460 

l M = 

9,2718 

L FAD I NG EDGE 

XL = 

-7.5514 

ZL = 

6 .5556 

LEADING SPIRAL COORDINATES 

XML = 

-13.8484 

ZMI_ = 

] 0.999 0 

TRAILING EDGE 

XM = 

5.5428 

ZM = 

8 .52 32 

trailing spiral coordinates — 

X M N= 

13.3810 

2 M M = 

1 1.7151 


X I ( i_ ) = 139.0378 XI I N ) = 56.3384 

* **V t R T I CA L REACT! ON *** 


WKPL= 4.6079 WKGI.= 1.8160 WKSL= -20.7413 TOTAL LEADING = -14.3173 

WKPT= 4.0207 WK.GT = 3.9375 WKST= -29. 5777 TOTAL TRAIL IN(>= -20.6194 

***HDK I 7 OUT A L REACTION*** 


HK PL = 17.0385 HKL = 12.3711 TOTAL. LEADING = 29.4096 

HKPT = -11 .7633 HK.T = -17.6463 TOTAL TrAIlING= -29.4096 


***Mt)MfcNTS**« 

MGL = -3.1632 MSI = 117.5277 MPL =-200.4954 TOTAL LEADING = -86.1314 
Mf,T = 18.9818 AST = -44.7006 MPT = 166.9795 TOTAL TRAILING= 141.2607 


TOTAL TORQUE = 55.1293 LB-1NCHFS = 7.6219 KG-MFTERS 

FM.FT, = .1968 0,0000 
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Fig. 3. Driven rigid roller 100% slip (s^. = 0, 
Vq = 0) on a pack of aluminum rolls 




Fig. 4. Soil-roller plastic flow configuration 
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Fig. 5. Soil-roller rim-velocity boundary conditions (bifurcation of plastic zones) 
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3 directions at leading and trailing zones 



14. Roller sinkage 
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